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ABSTRACT 

Spectra of galaxies contain an enormous amount of information about the rel- 
ative mixture of ages and metallicities of constituent stars. We present a compre- 
hensive study designed to extract the maximum information from spectra of data 
quality typical in large galaxy surveys. These techniques are not intended for de- 
tailed stellar population studies that use high quality spectra. We test techniques 
on a sample of globular clusters, which should consist of single stellar popula- 
tions and provide good test cases, using the Bruzual- Chariot 2003 high resolution 
simple stellar population synthesis models to simultaneously estimate the ages 
and metallicities of 101 globular clusters in M31 and the Magellanic Clouds by 
fitting their integrated spectra and photometry. The clusters cover a wide range 
of ages and metallicities, 4 Myr < tags < 20 Gyr and —1.6 < [Fe/H] < +0.3, es- 
timated by other methods in the literature. We compare results from model fits 
to both the spectra and photometry and find that fits to continuum-normalized 
spectra over the entire range available, typically 3500 A to 1 /im for this sample, 
provides the best results. For clusters older than 1 Gyr we agree with literature 
ages to 0.16 dex (35%) and [Fe/H] to 0.12 dex. For younger clusters we agree 
with literature ages to 0.3 dex (63%), but cannot constrain the metallicity. It is 
particularly important to use the entire continuum-normalized spectrum to avoid 
problems with model continua for young objects and to break age-met allicity de- 
generacies of broadband photometry. Our required S/N is 15-30 for 20% age 
uncertainties and 30-55 A^^ for 10% uncertainties over spectral resolutions of 
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AA = 5-25 A. This technique should work well for the age-metallicity parameter 
space expected for early-type galaxies at z~l, although individual galaxy spec- 
tral S/N may require the coaddition of a few like objects. Lack of accurate flux 
calibration in large surveys is not an issue for the continuum-normalized spectra. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — 
galaxies: individual (M31) — galaxies: star clusters — galaxies: stellar content 
— Magellanic Clouds 

1. Introduction 

Understanding when and how different types of galaxies formed in the universe is an age 
old problem. In recent years the hierarchical merging scenario (White & Rees 1978; Peebles 
1980; White & Prenk 1991), where large galaxies are built up over time through mergers of 
smaller objects, has been the leading theory for the mechanism by which structure formed in 
the universe. Large galaxy surveys have been undertaken to provide observational evidence 
for how galaxies formed. From their observations, wc are compiling a snapshot of galaxies 
througout time. At z~3 we find both massive starbursting galaxies in highly clustered 
regions (Neri et al. 2003; Blain et al. 2004; Chapman et al. 2004) and actively star forming 
galaxies elsewhere (Papovich, Dickinson, & Ferguson 2001). At z~2 we sec galaxies that are 
still forming stars, but have redder colors and larger masses from continued star formation 
since z~3 (Steidel et al. 2004). At z~l-2 a population of extremely red objects (EROs) are 
seen by near infrared surveys such as FIRES (Labbe et al. 2003), K20 (Cimatti et al. 2002), 
and MUNICS (Drory et al. 2001), many of which appear to contain older stellar populations 
that have experienced passive evolution (Forster Schreiber et al. 2004; Daddi et al. 2002; 
Saracco et al. 2003). At z~l we see a decrease in global star formation rates (Lilly et al. 
1996; Madau et al. 1996, 1998; Steidel et al. 1999; Barger et al. 2000; Thompson ct al. 2001; 
Wolf et al. 2005), an increase of stellar mass densities in early-types (Hogg et al. 2002; Bell 
et al. 2003; Dickinson et al. 2003; Kauffmann ct al. 2003b; Rudnick et al. 2003; Cimatti et 
al. 2004; Wolf et al. 2005), and a change in galaxy morphologies to more ordered systems 
(Strateva et al. 2001; Dickinson et al. 2003). At z<l we see early-type galaxies that appear 
to have formed at z~1.5-2 in the field (Im et al. 2002; Gebhardt et al. 2003) and at z~2-3 in 
clusters (van Dokkum et al. 1998; Thomas et al. 2005), possibly from those that were seen 
as star forming galaxies at z~2 and starbursts at z~2-3. 

This observed galaxy demography, however, is not matched by the semi-analytic hier- 
archical galaxy formation models. Although the models do well on predicting the global 
star formation rate and stellar mass densities, they do not match observations on the cor- 
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rect proportions of galaxy types as a function of redshift (Somerville 2004). In particular, 
not enough EROs or submm galaxies are formed at z~l-3 in the models. We do not yet 
understand the details of how galaxies assembled their mass, although the epoch of z~l-2 
continues to stand out as an important transitional period. Detailed studies of individual 
galaxies around this epoch will uncover essential clues about how they assembled into the 
distribution of galaxies that we see today. 

One such detail for study is determining ages or formation redshifts of galaxies in this 
epoch. Various techniques have been used in the past to determine an overall age from 
integrated galaxy light, including using colors to estimate age from the amounts of old 
(red) and young (blue) stars, using the equivalent widths (EWs) of specific spectral lines 
or breaks that are known to be sensitive to stellar age (Trager et al. 2000a,b; Kauffmann 
et al. 2003a,b), matching galaxy spectra to stellar spectra (Stockton et al. 1995; Dunlop et 
al. 1996; Spinrad et al. 1997), using fits of stellar population synthesis models to spectral 
energy distributions (SEDs) obtained from broadband or narrowband photometry (Papovich, 
Dickinson, & Ferguson 2001; Cimatti et al. 2004; Drory et al. 2004; McCarthy et al. 2004), 
and fitting galaxy spectra with model spectra (Yi et al. 2000; Nolan et al. 2003; Cimatti 
et al. 2004; Ferrcras & Yi 2004; McCarthy et al. 2004). We will refer to photometric SEDs 
throughout the paper as "SEDs," not to be confused with spectra. 

Photometry has been widely available and is a good way to cover a larger spectral 
range than typically available with spectroscopy to get information about different stellar 
populations, but it lacks spectral details and suffers from more degeneracy between age, 
metallicity, and dust content. Recent surveys with spectra of hundreds or thousands of 
galaxies, in addition to the photometry, provide more information for detailed analyses. 
Furthermore, the release of population synthesis models at high spectral resolution (Bruzual 
& Chariot 2003) allows one to take advantage of the detailed information contained in 
individual lines over the entire observed spectral range, rather than relying on only a few 
lines or on more degenerate broadband colors. 

In this work, we evaluate age and metallicity estimation techniques on globular clusters 
by comparing results from Bruzual & Chariot (2003) simple stellar population (SSP) model 
fits to the spectra, broadband photometry, and line indices of the clusters. Globular clusters 
provide a simpler test case than galaxies since their stars likely provide a coeval and nearly 
homogeneous metallicity population that has no internal dust extinction. We test the utility 
of using these models to estimate ages and metallicities of high redshift galaxies by first 
applying them to a sample of extragalactic globular clusters from the Santos et al. (2002) 
public database that have integrated spectra, photometry, and age or metallicity estimates 
obtained by other methods. Galaxies will be discused in a future paper. 
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This paper is arranged as follows. Globular cluster spectra and photometry are described 
in §2, comparison ages and metallicities from the literature are described in §3, the stellar 
population synthesis models are described in §4, our model fitting procedures are explained 
in §5, age and metahicity results are presented in §6, a discussion of issues is given in §7, 
and conclusions are summarized in 58. 



2. The Cluster Data 

The Santos et al. (2002) online spectral database provides our sample of extragalactic 
globular clusters that cover a significant range of age and metallicity, making their results 
applicable to evolving galaxies. We chose those clusters that also have photometry available 
in the literature, which includes 79 clusters in the Large Magellanic Cloud (LMC), 4 clusters 
in the Small Magellanic Cloud (SMC), and 18 clusters in M31. The spectra were obtained 
with the 1.52 and 2.2-meter telescopes at ESQ in La Silla (Bica et al. 1986a, 1987a,b, 1990, 
1994), the 2.15-m CASLEO telescope in Argentina (Santos et al. 1995), the 3.6-m CFHT 
in Hawaii (Jablonka et al. 1992), and the 4.2-m WiUiam Herschel Telescope in La Palma 
(Jablonka et al. 1998). 

The data include near ultraviolet, optical, and near infrared spectra, covering a total 
wavelength range of 3200-10000 A for some clusters, and varying in spectral resolution from 
6 to 23 A. We should note that for all clusters except a few in M31, this resolution element is 
larger than the AA = 10 A required for Lick/IDS line indices. All available spectral segments 
for each object are used to utilize maximum wavelength coverage, but since they originate 
from a number of different sources, flux calibration between regions may not be consistent. 
Table 1 lists each object with its specific wavelength coverage, spectral resolution, and source 
references. 

Wc compile photometry of the clusters from various sources in the literature. Bica et al. 
(1996) provide UBV photometry for the LMC clusters with ~20% of the objects taken from 
van den Bergh (1981) and the rest rc-obscrvcd by the authors at CTIO or CASLEO. Persson 
et al. (1983) provide VJHK photometry for LMC clusters from observations on 1-2.5 meter 
telescopes at Las Campanas and Cerro Tololo observatories. For the SMC clusters, Persson 
et al. (1983) provide VJHK and UBV are from Santos et al. (1995) who use photometry 
from van den Bergh (1981) and a few other sources. For the M31 clusters, Battistini et al. 
(1993) supply BVRI from observations on the 152-cm telescope in Loiano, Bologna and on 
the 4-m KPNO telescope. Barmby et al. (2000) provide UBJHK from observations on the 
1.2-m telescope at the Fred L. Whipple Observatory. 
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3. Compcirison Ages and Metallicities 

Even for a given globular cluster, a wide range of age and metallicity estimates typi- 
cally exists in the literature. This is not completely surprising given the different methods 
employed to derive the parameters, some based on broadband colors and others on very 
narrowband spectral line strengths, but the dispersion in estimates is quite large for many 
clusters. We will compare our results to age and metallicity values from a number of sources 
and different techniques to show where our values lie in the total spread. Literature ages 
and metallicities for specific clusters can be found in Table 4, along with our derived values. 

One of the oldest techniques for measuring the age of a star cluster is using its integrated 
UBV colors. When placed on color-color plots, globular clusters follow a sequence that can 
be related to age. We use MC cluster U-B color ages that were determined by Bica et al. 
(1990) and Santos et al. (1995). When photometry is available for individual stars in a 
cluster, the location of the main sequence turnoff (MSTO) on a color-magnitude diagram 
(CMD) can be used to determine the age of the cluster. We use a large compilation of MSTO 
ages for MC clusters from Hodge (1983). Refinement of the MSTO technique is achieved by 
fitting theoretical isochrones to the CMDs. We take CMD fit ages from Elson & Fall (1988) 
and Girardi et al. (1995, 1998). We also add comparison ages determined from photometric 
fits to GISSEL96 model SEDs (Bruzual & Chariot 1993) in the BATC filter system, using 
15 intermediate bandwidth passbands, by Jiang et al. (2003), and from the line index ratios 
R5 / Fel A4045 and Ca II H/K by Leonardi & Rose (2003), which were developed to work 
better on younger objects than the typical Lick line indices for use on older objects. We 
consider the Leonardi & Rose line index ratio ages to be some of the most secure, given 
that they are determined from high quality spectra and are tailored for the younger ages of 
many clusters in this sample. Santos & Piatti (2004) utilize diagnostic diagrams constructed 
from EW sums of metal (Ca II K, G band, Mg I) and Balmer (US, H7, H/?) lines, which are 
calibrated against literature values and placed on homogeneous age and metallicity scales, 
to determine ages and metaUicites of MC clusters. We also include these in our comparison 
values of age and metallicity. 

Fitting theoretical ishchrones to cluster CMDs is generally accepted as the most secure 
age determination possible when using photometry, however the results do vary between sets 
of ishochrones. The absolute derived age depends on model zeropoints and uncertainties. 
Schiavon et al. (2002) find that Padova isochrones with AGB stars included underpredict 
the luminosity function of giants by 0.2-0.3 dex when compared to observations of 47 Tuc, 
resulting in an overprediction in the spectroscopic age of 2-3 Gyr for old stellar populations. 
Young stellar populations are much less affected. Additionally, the Padova isochrone uncer- 
tainty in Teff of 75 K leads to a 1 Gyr age uncertainty or 0.1 dex uncertainty in the [Fe/H] 
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scale. Vazdekis et al. (2001) also find discrepancies between CMD and spectroscopic ages of 
47 Tuc due to the exclusion of a-enhancement and atomic diffusion in evolutionary models, 
causing a zeropoint offset. If 47 Tuc is representative of other old stellar populations, for 
clusters in the 10-14 Gyr age range these model isochrone issues could lead to spectroscopic 
based age discrepancies of 20-40 % when compared to CMD based values. 

Large homogeneous samples of globular cluster mctallicitics arc harder to find, since 
metallicity cannot be as accurately determined from photometry and typical methods require 
high quality spectra. However, some photometric techniques have been calibrated to spectra. 
Of these we use [Fc/H] values from VJK colors by Cohen & Matthews (1994), fits to CMDs 
(Seggewiss k Richtler 1989; Sagar k Pandey 1989; Suntzeff et al. 1992; Girardi et al. 1995; 
Dirsch et al. 2000; Johnson et al. 2001; Piatti et al. 2002), and fits to GISSEL96 model SEDs 
by Jiang et al. (2003). Our more secure comparison metallicites come from spectra-based 
techniques. For these we use [Fe/H] from the Leonardi & Rose (2003) fine index ratios and 
their comparison literature values that were derived from EWs of Ca, Fe, Mg, and Na lines 
(Cohen 1982), and from EWs of the Ca II triplet at A8500 (Olszewski et al. 1991; Da Costa 
k Hatzidimitriou 1998). We incorporate the comparison literature values used by Santos k 
Piatti (2004), which include metal line abundances (Jasniewicz k Thevenin 1994; Hill et al. 
2000), line indices (Beasley et al. 2002), and an infrared index at 1.6 /xm (Oliva k Origlia 
1998). And finally we use a combination of line indices and spectral breaks from Huchra et 
al. (1991): D(4000), CNB A3883 A, G-band A4300 A, MgH, Mg b, and Fe A5270 A. 



4. Stelleir Population Synthesis Models 

We use the high spectral resolution (R~2000) SSP "standard model" from Bruzual k 
Chariot (2003; hereafter BC03) that utilizes the STELIB/BaSeL 3.1 spectral library, the 
Padova 1994 evolutionary tracks, and the Chabrier (2003) initial mass function (IMF) with 
stellar mass fimits of 0.1 and 100 Mq. They provide model spectra at [Fe/H] = —2.25, —1.65, 
—0.64, —0.33, -1-0.093, and -1-0.56, which we finearly interpolate at each wavelength point to 
a grid of spectra that spans [Fe/H] of —2.2 to -1-0.5 in increments of 0.1 dex, and ages of 1 
Myr to 20 Gyr in increments of 2 Myr for 1< tage <10 Myr, 20 Myr for 10< tage <100 Myr, 
200 Myr for 100< tage <1000 Myr, and 1 Gyr for 1< tage <20 Gyr. 

The BC03 models include thermally- pulsing asymptotic branch stars (TP-AGB), mak- 
ing use of the multi- metallicity models of Vassiliadis k Wood (1993) that have been calibrated 
on stars in the Galaxy, LMC, and SMC. The importance of these stars, which have a strong 
influence on the integrated near infrared light from star clusters of certain ages, will become 
apparent later in the paper. This phase can dredge up carbon in the stellar atmospheres. 
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leading to carbon-rich stars, which are also very red and can dominate the near infrared hght 
from some clusters. Although no simple prescription can be expected to match all clusters, 
these stars must somehow be included. Bruzual & Chariot (2003) achieve this by defining 
the transition to carbon stars and the duration of this phase with the models of Groenewe- 
gen & dc Jong (1993) and Groenewegen et al. (1995). This semi-empirical prescription for 
TP-AGB and carbon stars has been tested on and provides good agreement with observed 
colors of Magellanic Cloud clusters and with optical and near infrared surface brightness 
fluctuations of metal-poor Galactic globular clusters and more metal-rich nearby elliptical 
galaxies (Liu et al. 2000). 

5. Data Prepairation and Model-Fitting Procedures 

Our approach is to use the maximum amount of information possible for each object. 
Broadband photometry often covers a much wider spectral range than spectroscopy, provid- 
ing better constraints on the possible mixes of stellar types, but suffers from age-metallicity 
degeneracy. However, by combining broadband information with detailed spectral analysis, 
we hope to better constrain model flts and thus derived ages and metallicities. We investi- 
gate the ability of both broadband and detailed spectral features of the models to estimate 
age and mctallicity of the globular clusters by fitting models to the full spectrum, to the 
continuum-normalized spectrum, to the continuum shape, to photometry, and to spectral 
line indices of the clusters. Each of these fitting procedures is described in detail in the 
following sections. 

• Full spectrum fits utilize information from both the line strengths and the continuum 
shape. 

• Continuum-normalized ( CN) spectrum fits use only the information contained in the 
lines, losing important continuum information, but also removing any adverse effects of 
inaccurate ffux calibration in the data or continuum shape errors in the models. Model 
continua problems are particularly possible in the near infrared due to the inability to 
properly account for the number of TP-AGB or carbon stars that can contribute much 
of the light at these wavelengths. 

• Continuum fits use only the continuum shape that was removed from the spectrum. 

• Photometry fits cover a much broader wavelength range, but lose detailed spectral 
information and may also be affected by continuum shape problems in the data or 
models. 
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• Line index fits move from broad to narrow and focus only on specific spectral lines 
that are known to be sensitive to metallicity or age. We evaluate these specific lines 
in the models by calculating the indices directly from model spectra and comparing 
those to line indices calculated from the cluster spectra in the same way. 

5.1. Spectra 

The spectra for each object in the Santos et al. (2002) database consist of near ultra- 
violet, optical, and near infrared segments from multiple sources at spectral resolutions of 
6-23 A. We splice all spectra together to cover the widest possible range for each cluster. 
No noise or sky spectra are included in the database, so we must estimate the noise from 
the spectra. To do this we calculate the standard deviation of the spectrum in 100 A bins 
and scale the resultant signal-to- noise ratio (S/N) spectrum to match the average S/N value 
quoted in the source publications for some objects. For those objects that do not have quoted 
S/N values, the standard deviation spectrum provides the noise. This hkely overestimates 
the noise since absorption lines will increase the standard deviation of a bin, but the required 
scahng is typically less than 10%. 

We prepare the data for fitting by marking regions that contain emission lines, noise 
spikes, and sky background residuals, as well as two regions for which the STELIB stellar 
hbrary used in the models has problems with telluric features (6850-6950 A and 7550-7725 A). 
These regions are masked and deweighted for each cluster by significantly lowering its S/N 
spectrum over the affected wavelengths. We smooth the grid of model spectra, which have an 
intrinsic resolution of 3 A, to the wavelength dependent resolution of each cluster and fit its 
spectrum using two free parameters, age and metallicity, and a S /N-weightcd normalization 
constant, an, which is uniquely determined for each model. We select the best fitting model 
by calculating the summed for estch model spectrum of a given age and metallicity, relative 
to the cluster spectrum, and selecting the model in the grid with the minimum value. In 
the calculation, we use the previously determined cluster noise spectrum as the (7j for the 
flux, fi, at each wavelength bin, i, for each model, n, as given in Equations 1 and 2. Once 
the best fitting model is determined for a cluster, wc inflate its noise until the xleduced—^ 
to include any additional noise sources and to better estimate confldence contours and error 
bars. If, after this process, the estimated errors are smaller than the distance to the next 
model grid point, an error bar of half this distance is assigned to the cluster for both age 
and metallicity. 
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Errors on age and metallicity estimates are determined from the A^^ contours. The 
S/N of the data controls the size of these contours, and thus the confidence of the derived 
parameters. Typical levels for the globular clusters are S/N~100 per resolution element. 
Figure 1 shows a simulation of the effect of S/N on the confidence of the derived age and 
metaUicity for NGC 419 in the SMC, achieved by artificially adding different levels of gaussian 
noise to its spectrum and fitting models to each of the noisy continuum-normalized spectra. 
The simulated S/N per resolution element from top left to lower right in Figure 1 is 99, 
22, 12, 8, 6, and 4. Solid, dashed, and dotted lines represent 1, 2, and 3cr contours. The 
color bars give the Ax^ values for confidence levels of 68.3% (Ax^=2.3), 90%, 95.4%, 99%, 
99.73%, and 99.99% (Ax^=18.4) for two degrees of freedom. Black dots mark the locations 
of models in the age-metallicity grid. The best fitting model is marked by a black cross in 
each panel. The age and metaUicity from the best fitting model in the highest S/N case is 
marked by a yellow circle in each subsequent panel for comparison. 

Two things become apparent from these plots. First, as the confidence contours grow 
with decreasing S/N, metallicity becomes harder to constrain than age. In some low S/N con- 
ditions (e.g. the last panel) multiple islands of age-metallicity combinations arise as equally 
likely answers. The best fitting model in each panel does not always overlap the la con- 
tours of all other panels, however, this mostly happens in the [Fe/H] dimension and those 
values are less secure than the ages. This is our first indication that metallicity is poorly 
constrained, compared to ages. The second point is that we find a limiting S/N below which 
the estimated age and metallicity become unreliable. This happens somewhere between S/N 
= 22 and 12 in these plots, at which point the contours in the low S/N cases do not overlap 
the estimates from the highest S/N case. We find that the S/N must be >15 per resolu- 
tion element (S/N~5 A~^) for AA=13-23 A to adequately constrain the age. MetaUicity 
is harder to constrain under all conditions. A more detailed analysis of limiting S/N and 
spectral resolution can be found in §6.5. 
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5.2. Continuum-Normalization 

Errors in flux calibration of the spectra can induce errors in derived ages and metal- 
licities. Additionally, model continua can be off in the near infrared due to stochastic con- 
tributions from TP-AGB stars in ~ 0.1 to 1 Gyr old populations and from carbon stars in 
populations of age 0.3 to 2.5 Gyr (Frogel et al. 1990; Marigo et al. 1996; Girardi et al. 1998). 
These variations can result in differences of nearly 2 magnitudes in V-K (Bruzual & Chariot 
2003, Figure 8). AGE stars in the thermally pulsing phase can contribute over 80% of the 
K-band light (Maraston 2005), the exact amount of which is hard to predict since this phase 
involves ejection of the outer stellar envelopes and geometry-dependent obscuration of the 
central stars. The BC03 models include TP-AGB and carbon stars semi-empirically, but in 
any given cluster the K-band light can vary depending on the numbers of stars it actually 
has in these phases. 

To avoid these issues, we normalize out the continua with a median-binning routine 
that uses bin sizes tailored to work for different classes of objects. The young clusters of the 
Magellanic Clouds have strong Balmer discontinuities and need to be sampled by small bins 
to correctly trace this sharp continuum break, while older objects in M31 require larger bins 
to smooth over absorption troughs in the red part of the spectrum. The bins that we use are 
100 A wide blueward of 4000 A and 300-700 A wide redward of 4000 A. The points in each 
bin are sorted by flux values, the lower 1/2 of the points are ignored to reduce the weight of 
absorption lines, the highest few points are rejected to ignore noise spikes, and the median 
of the remaining flux values is assigned to that bin. The medianed bin fluxes are linearly 
connected and the object spectrum is divided by the result, providing a flat continuum. The 
same binning parameters used for an object are also used on each model spectrum in the 
grid. 

Removal of the continuum shape should not affect the values of line indices measured on 
the spectrum, providing a good check of our procedure. We applied this test, as suggested 
by the referee, by choosing a random sample of 25 globular clusters spanning the entire age 
range and measured line indices before and after continuum normalization. The indices used 
were Mg6, Fe5270, Fe5335, H/9, RSa, H(5f, H-fA, ^If, G4300, <Fe>, and [MgFe]'. Only the 
positive index values that could be reliably measured were used, as in §6.4. On average the 
values change by less than 1%, and all line index median differences for the 25 clusters are 
below 3% with no correlation to the derived cluster age. 

The continuum shape that we divide out is also fit by model continua and referred to 
throughout the paper as "continuum fits." Different types of spectral fits are illustrated in 
Figure 2 for an old cluster, G177 in M31, and a young cluster, NGC 1711 in the LMC. 
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5.3. Photometry 

To test ages and metallicities derived from broadband features over a wide baseline, we 
fit models to the UBVJHK photometry from the literature, using all of the bands that are 
available for each cluster. Our model fitting procedure requires scalable fiuxes, so we convert 
cluster magnitudes to fluxes using the UBV zeropoints from Bessell (1979) and the JHK 
zeropoints from Wamsteker (1981). We obtain model fluxes for each band by convolving 
fllter transmission curves with the model spectra. These fluxes are scaled to those of the 
object with an overall S/N-weightcd normalization parameter during the fit. The best fitting 
model is selected by minimization, using Equations 1 and 2 with i now representing 
each photometric band. Errors are not given for all literature photometry, so we initially 
assume a 5% photometric fiux error in all bands for each cluster, which is then adjusted 
to make Xreduced~^ the best-fitting model. Fitting parameters are age and metallicity. 
Photometric fits are also illustrated in Figure 2. 

5.4. Line Indices 

Lick line indices have been used for some time (Faber et al. 1985; Burstein et al. 1986; 
Gorgas et al. 1993; Worthey et al. 1992, 1994; Worthey 1994; Worthey & Ottaviani 1997; 
Trager et al. 1998). They are based on a very specific set of data and their application to 
new objects is somewhat complicated. To properly use Lick/IDS indices, some of the same 
stars must be observed with the instrumental setup of the new objects and the spectra of the 
new objects must be convolved to exactly the same wavelength dependent resolution as the 
original Lick data, which varied from run to run during the development of the indices. Very 
small differences in the wavelengths can cause large errors in the indices. To investigate the 
robustness of these specific lines in the model spectra, we calculate line indices directly from 
the models, similar to the calibrations done in Bruzual & Chariot (2003). Direct application 
of the models for this purpose would avoid the step of convolving the object spectra to the 
estimated Lick resolutions. If the models are smoothed to the same spectral resolution as 
the cluster data, indices can be calculated directly from both spectra in the same manner. 
However, the bandpasses of the Lick indices were originally chosen for data with AA = 8- 
12 A (Worthey & Ottaviani 1997). Specific lines in lower resolution spectra may begin to 
lose age and metallicity information. Nevertheless, given the long history of the use of line 
indices, we wanted to attempt a variation of this technique on data with resolution and S/N 
representative of that from typical galaxy surveys. 

First, we broaden the model spectra to match the resolution of each object, which may 
vary for wavelength regions obtained on different instruments. Then, we calculate indices 
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for the lines in the Lick system using the passbands defined in Trager et al. (1998) and 
Worthey & Ottaviani (1997). To these we add <Fe>, an average of Fc A5270 and Fe A5335; 
D„(4000), the 4000 A break strength using the narrow passbands defined in Balogh et al. 
(1999); CNB A3883 A, MgH, and Ca II H+K, as defined by Brodie & Huchra (1990); and 
[MgFe]' = v/Mg6(0.72 x Fe5270 + 0.28 x Fe5335), as defined by Thomas ct al. (2003). The 
[MgFe]' index has been found to be a good tracer of total metallicity and to be insensitive 
to cu/Fe enhancement (for spectral resolution near 10 A), which may be important for some 
of the globular clusters in this sample. 

We calculate the indices as described in Trager ct al. (1998), by connecting the average 
flux values of the pseudocontinua sidebands with a straight line and integrating the fiux 
under that line over the index passband. Indices (or break strengths) are calculated for 
CNB, RSa, HSf, CNi, CN2, Ca4227, G4300, H7A, H7F, Fe4383, Ca4455, Fe4531, C24668, 
H/3, MgH, Fe5015, Mgi, Mg2, Mg&, Fe5270, Fe5335, Fe5406, Fe5709, Fe5782, Na D, TiOi, 
Ti02, [MgFe]', <Fe>, and D„(4000). Indices for model spectra are calculated in the same 
manner. We select the model with the smallest sum of absolute value residuals between 
its index values and the corresponding indices derived from a cluster's data. The age and 
metallicity of this closest matching model are adopted as the estimates for the cluster; no 
interpolation is performed between model grid points. To facilitate comparisons with the 
hterature, we repeated the model selection using different subsets of the indices hsted above: 
all indices simultaneously; Huchra et al.'s (1991) group of 6, CNB, D„(4000), G4300, MgH, 
Mg&, and Fe5270; and H/3-[MgFe]' as in Thomas et al. (2003). 



6. Results 

Examples of model fits are shown in Figure 2 for old and young clusters. Panels a,c,&e 
are the cluster G177 in M31 and panels b,d,&f are NGC 1711 in the LMC. The fuU spectrum 
fits are shown in a&b, the continuum-normalized spectrum fits are shown in c&d, and the 
continuum fits are shown in e&f. Photometry fits are also shown in a&b with U, B, V, J, 
and K bands for G177 and U, B, V, J, H, and K bands for NGC 1711. The circles are the 
object photometry and triangles are model fiuxes from the best photometric fit. The solid 
lined black spectra are the clusters, the dashed cyan lines are best fitting models from the 
spectral fits, and the dotted green lines are spectra corresponding to the best fitting models 
from the photometry fits. The model spectra corresponding to the best fitting photometry 
are normalized to the cluster spectra at V band in these plots. In the continuum fits, black 
lines are the data and green lines are the models. The resulting ages and metallicities from 
the different types of fits are given in the plot labels. 
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For these two clusters the ages derived from the different types of model fits nearly agree 
at 20 Gyr for G177 and about 0.1 Gyr for NGC 1711. The agreement between photometric 
and spectroscopic fits seen here is not always the case. G177 in Figure 2a illustrates a case 
where the photometry through K band reveals a redder SED than the optical spectrum. In 
this case, the best fitting photometry model has a higher metallicity than the model that 
best fits the spectrum to make the object redder, an effect of age-metallicity degeneracy. 
A similar effect is seen in some of the young clusters, presumably due to a near infrared 
excess from TP-AGB or carbon stars in that case. Some examples of such clusters within 
the affected age range of 0.1-2 Gyr are shown in Figure 3. 



6.1. Different Model Fitting Techniques 

To evaluate the different model fitting techniques, we first compare their results to each 
other. In the following sections we compare our results to the literature. Figure 4 compares 
the ages derived from CN spectrum, full spectrum, photometry, and continuum fits. Clusters 
with average literature ages > 1 Gyr are maked with yellow circles and < 1 Gyr with cyan 
diamonds. Although the general trends of the ages agree, there are differences in the results. 
When compared to CN spectrum fits, the methods that contain the continuum overestimate 
the ages of many clusters younger than 1 Gyr. Wc believe that this is due to red supergiant 
stars for clusters in the age range of a few to tens of Myrs and TP-AGB stars for clusters 
in the age range of 0.1 to 1 Gyr. Both of these types of stars can have stocastic effects in 
star clusters because of their extremely high luminosities. If more light from them appears 
in the integrated spectra of the cluster than is included in the models, the cluster SED will 
appear redder than the model, forcing an older redder model SED as the best fit. All fitting 
methods that contain model continua will have this same bias. This effect is discussed further 
in comparison with literature ages in §6.2. 

For clusters with average literature ages < 1 Gyr our tightest estimated age correlation 
is seen between the full spectrum and continuum fits (Figure 4c), both of which contain the 
spectral shape. The median fractional age offset, relative to the average literature age, be- 
tween these fits is 0.04 for 63 clusters. The fact that the continuum-full spectrum correlation 
is tighter than the CN-full spectra correlation (fractional offset of 2.40) suggests that the 
continuum shape has a larger effect than the spectral lines in constraining the ages of these 
young globular clusters for our employed fitting methods. Therefore, any problems in the 
continuum shape of the data or models will have grave effects on the derived ages. For the 
clusters older than 1 Gyr our tightest estimated age correlation is between the CN spectrum 
and photometry fits (Figure 4b), with a median fractional offset of 0.18 for 12 clusters. 
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For young clusters the metallicities obtained from our different techniques show virtu- 
ally no correlation, with a tendency for the continuum to give higher metallicities than the 
full spectrum. For clusters older than 1 Gyr there is a very weak correlation between metal- 
licities obtained with continuum and full spectrum fits, possibly indicating that the derived 
metallicity has a stronger dependence on the continuum shape than the spectral lines for 
our fitting method. 

6.2. Ages 

In Figure 5 we compare our globular cluster ages to those obtained by other methods in 
the literature. The thick cyan horizontal bars connect points that represent the same cluster, 
but that have ages from multiple sources in the literature. These illustrate the large spread 
in previously derived ages. The thin horizontal bars are the quoted errors on the literature 
ages and the thin vertical bars are our la errors on ages derived from model fits. Note that 
the errors are much larger for the photometry fits, reflecting more degeneracy when using 
broadband colors than when using spectra or the continua shapes. 

Over the entire age range spanned by the clusters and considering all literature age 
estimation techniques, our CN spectrum fit ages best match those from the literature. The 
average offsets and rms scatter our ages from groups of literature values based on different 
techniques are given in Table 2, where column 1 gives values for all clusters in the sample that 
have literature ages and column 2 excludes two outliers (for reasons explained in §6.2.1) and 
the clusters that only have literature values from photometric fits to BC96 models (because 
these used only a few fixed metallicity values in the model grid, which could skew the 
derived ages, and they are based on old models). Excluding these points, our CN spectrum 
fits compared to the averages of the literature values for each object have an average offset 
of 0.29 dex (0.19 Gyr) with a dispersion of 0.36 dex (1.00 Gyr). 

Our photometric ages have an average offset of 0.59 dex (0.48 Gyr) and dispersion of 
0.76 dex (1.19 Gyr). There arc some systematic differences between the two sets of results. 
The scatter of the CN spectrum fit points is more uniformly distributed around the literature 
values, unlike the photometric points that tend to overestimate the age, particularly for the 
younger clusters. The full spectrum fits, which contain the continuum, also overestimate the 
ages below 1 Gyr, with an average offset of 0.69 dex (0.54 Gyr) and a dispersion of 0.9 dex 
(3.2 Gyr). This is the same effect that was seen when comparing our fits with and without 
the continua included. 

We now return to stars that can strongly affect the continua shapes. Detailed studies 
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of such stochastic effects in star clusters can be found in Bruzual & Chariot (2003), Cervino 
et al. (2000, 2001, 2002, 2006), and Maraston (2005). Our photometric fits, in particular, 
show two different age ranges where our best-fit ages depart from literature ages on the high 
side (Figure 5c): 4 to 100 Myr and 0.1 to 1 Gyr. The age range of 0.1 to 1 Gyr would be 
affected by TP-AGB stars. If the model prescription for these stars does not exactly match 
those in the clusters, then the older ages from fitting methods that include the continuum 
could be explained by older, redder spectra giving the lowest for clusters that have this 
near infrared excess. Our derived ages also scatter to older values for clusters younger than 
0.1 Gyr when using models that include the continuum shape. This effect is hkely due to 
stocastic red supergiants, which have been observed with estimated masses up to 120 Mq in 
the LMC (Massey et al. 2005). These stars will affect clusters younger than a few tens of 
Myrs. Either the number of red supergiants is underestimated in the BC03 models, or the 
upper mass cutoff of 100 Mq for the Chabrier IMF is not quite high enough for the LMC. 
See §6.2.1 for a discussion of two LMC clusters that are dominated by supergiants, which 
are also our extreme outliers in the CN spectrum fits of Figure 5a. 

Because of the large age range covered by the clusters, it is perhaps useful to look at 
the behavior of the fractional age errors relative to the literature values, which are plotted 
in Figure 6. Overall, the CN spectrum fits provide ages with fractional offsets from the 
hterature of 61%, while the photometry does much worse at 1048%. If we consider younger 
and older clusters separately, the CN spectrum fits do better than photometry on the average 
age offsets for both young and old clusters, but the photometry does better on fractional 
errors for the older objects. For clusters with ages < 1 Gyr, the CN spectrum fits have an 
average offset of 63% (0.30 dex or 0.09 Gyr) and dispersion of 0.37 dex (0.35 Gyr), while 
the older clusters have an offset of 35% (0.16 dex or 1.63 Gyr) and dispersion of 0.17 dex 
(3.78 Gyr). The dispersion is higher below ages of 0.1 Gyr even though the CN spectra are 
not affected by the continua effects mentioned earlier. This increased dispersion is likely 
due to the degeneracy of many spectral lines in this age range (see Figure 10 and §6.4). 
The photometry fits do much worse fractionally on young clusters with an average offset 
of 1115% (0.62 dex or 0.38 Gyr) and dispersion of 0.78 dex (1.02 Gyr), and better on the 
fractional errors of older clusters with an offset of 23% (0.10 dex or 2.06 Gyr) and dispersion 
of 0.11 dex (2.67 Gyr). 

The lowest fractional errors for the CN spectrum fits occur for the Leonardi & Rose 
(2003) line index ratios and for the clusters older than 1 Gyr, both at 35%. Photometry fits 
do shghtly better on the older objects with errors of 23%, but these errors rise substantially 

for the young objects to 1115%. The largest fractional age offsets occur in the clusters less 
than 1 Gyr old for the full spectrum (3450%), continuum (4432%), and photometry fits. This 
is likely due to AGB stars. Furthermore, the full spectrum and continuum fits have larger 
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age offsets in the older clusters (425% and 480%) than do the CN spectrum or photometry 
fits, which may be due to flux calibration issues, since both of these methods include the 
observed spectral continuum shape. To take care of both data and model continua issues, it 
appears that continuum-normalized spectral fits are the most robust in producing accurate 
age estimates for the clusters of all ages. 

The age errors do not correlate with mctallicity, so they seem to be dominated by the 
age regime of the cluster. Although photometry does a decent job on older objects, it should 
be noted that after excluding the BC96 photometry literature points, we only have 4 clusters 
in the > 1 Gyr bin at ages of 1.2, 2.2, 12.4, and 13.1 Gyrs. It could just be that the two 
younger ones have no problems with boosted near infrared emission. In general, since there 
is not a good way to determine whether the 1-2 Gyr old objects are affected by the AGB 
star problem, using the CN spectrum fits would avoid this issue. Some combination of 
CN spectrum and photometry fits will likely provide the most robust answer. This will be 
explored further in future work. 

Other studies have compared results from their cluster age estimation techniques to 
the hterature. Rafelski & Zaritsky (2004) use colors formed from UBVI photometry along 
with Starburst99 (Leitherer et al. 1999) and GALEV (Anders & Pritze-v. Alvensleben 2003) 
models to derive ages and compare to those in the literature from integrated colors (van 
den Bergh 1981; Hunter ct al. 2003) and isochronc fitting (Pictrzyhski & Udalski 1999; de 
Olivcira ct al. 2000; Mighell ct al. 1998; Rich et al. 2000) in Figure 8 of their paper. They find 
an age correlation with dispersion of 0.76 Gyr, which drops to 0.49 Gyr when considering 
only the more secure literature ages derived from GMDs. Over a similar age range our CN 
spectrum fits provide ages with an overall dispersion of 1.0 Gyr about literature comparison 
values, with 0.12 Gyr for UBV colors and 1.06 Gyr for CMD fits. Our photometry fits have 
dispersions of 1.04 Gyr for UBV colors and 0.85 Gyr for CMD fits. Although Rafelski & 
Zaritsky (2004) use clusters in the SMC, only two overlap with our sample. Our hterature 
ages come from different sources as well, so we cannot make a direct comparison to this 
work. 

To summarize our results for globular cluster age estimation, the CN spectrum model fits 
best match the entire range of age estimates from the literature (Figure 5 and Table 2), with 
average errors of 0.16 dex (35%) for older clusters and 0.3 dex (63%) for younger clusters. 
The full spectrum, continuum, and photometry fits overestimate the ages of many clusters 
below 1 Gyr, apparently suffering from an excess of AGB stars making the cluster spectra 
redder than the models and forcing older aged models to provide the best fits. Although 
photmetric fits seem to do better for clusters older than 1 Gyr (0.1 dex or 23% error), we will 
show in §6.3 that CN spectrum fits are superior for simultaneously providing both age and 
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metallicity for these older objects. Furthermore, the uncertainty of whether the clusters will 
have boosted near infrared emission due to AGB stars makes using the photometry alone 
a less reliable technique than using spectra. Comparing the fitting methods (Figure 4), we 
find a tighter correlation between ages derived from the full spectrum and continuum fits 
than from the CN and full spectrum fits, suggesting that the ages are more strongly driven 
by the continuum shape than by the spectral lines. Because the derived ages are so strongly 
influenced by the continua and we have seen signs of problems in the models matching the 
cluster continua, we conclude that the best method for deriving accurate ages for globular 
clusters of all ages, especially when simultaneously determining metallicity, is fitting models 
to their continuum-normalized spectra. 



6.2.1. The Outliers 

There are two extreme outliers in the CN spectrum age plot of Figure 5a. The spectra of 
these outlying objects are shown in Figure 7. The upper leftmost outher in the CN spectrum 
age plot is NGC 2092 in the LMC (Figure 7a,b), with a literature age of 4-12 Myr (from 
MSTOs and UBV colors) and for which we derive an age of 10 Gyr from both the CN and full 
spectrum fits. Its spectral coverage is only over the optical range of 3500-5870 A. The shape of 
its spectrum is fairly fiat with emission lines and few absorption features. The emission lines 
include [011]A3727, [Olll]AA4959,5007, and [Nelll]A3869 nebular emission. The emission 
lines are masked out in our fits, which leaves very few distinguishing features. The cluster's 
UBV photometry is also essentially flat. Our best flt metallicity of this object is bottomed 
out at the minimum in the model grid, [Fe/H]=— 2.2, for both the CN and full spectrum 
fits, likely because nothing would fit well and the lowest metallicity model that has very 
few metal fines in the red produced the lowest, although high, -y^. Santos et al. (1995) list 
this object as a cluster embedded in a star forming complex with a fiat continuum probably 
caused by red supergiants, which arc stocastic in nature even for large star clusters. This 
appears to be a case that cannot be fit well with either spectral or photometric techniques 
using the BC03 SSP models. 

The other outlier in the upper lefthand corner of the CN spectrum age plot in Figure 5a 
is NGC 2096 in the LMC. This object has a literature age of 49 Myr from its U-B color. 
Its spectrum (Figure 7c, d) contains a moderate sized 4000A break of D„(4000)=1.28, which 
is consistent with an age of at least 1 Gyr from index plots in Bruzual & Chariot (2003), 
and shows no evidence of the characteristic spectral shape of a young object with a strong 
Balmer discontinuity that would be expected for an age of only 49 Myr. For this object, our 
age is 1.2 Gyr (full spectrum and photometry fits) to 6 Gyr (CN spectrum fit). However, 
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Santos et al. (1995) note this object as one where a few luminous intermediate temperature 
supergiants dominate the cluster's integrated spectrum and claim that its age is even younger 
at 6-12 Myr. This may be a stocastic case where the supergiant dominated spectrum is not 
properly matched by the models. 

Two out of 101 globular clusters in this sample could not be properly fit by the BC03 

models. Both are very young clusters that show evidence of supergiant dominance in their 
integrated spectra. The stocastic nature of these stars make it difficult for us to derive 
accurate ages for clusters younger than a few tens of Myrs. 

6.3. Metallicites 

We compare our derived cluster metallicites to the literature values in Figure 8. There 
arc fewer points on these plots than for the ages in Figure 5 because fewer of the clusters 
had mctallicity estimates in the literature. The magnitudes of the metallicity errors relative 
to literature values are a stronger function of cluster age than [Fc/H]. This is illustrated in 
Figure 9 where the [Fe/H] errors arc plotted against average literature ages. Table 3 gives 
the average [Fe/H] offsets and dispersions for different groups of clusters. 

Metallicity is harder to constrain than age. This was first seen in Figure 1 when the 
Ax^ contours grew primarily in the metallicity dimension as noise was added to the cluster 
spectrum. Nevertheless, our CN spectrum fits do well on metallicity estimates for the older 
clusters. In Figure 8 we see a tight correlation of metaUicity from these fits with the hterature 
for clusters older than 1 Gyr. The large symbols mark literature estimates that are based on 
spectra. There is an overall [Fe/H] scale offset between M31 {blue pluses) and the MC {red 
circles, black crosses, black inverted triangles) . If this offset of 0.45 dex is removed from the 
MC cluster literature metallicities, our metallicites agree to 0.12 dex with those from the 
Huchra ct al. (1991) line indices {blue pluses), literature values on a homogeneous scale from 
Santos & Piatti (2004) {black inverted triangles), and Lconardi & Rose (2003) line index 
ratios {red circles) and their comparison literature values {black crosses). 

The small symbols have literature metallicity estimates that are based on colors. There 
is more scatter in these values. The weaker correlation with VJK colors {green triangles) 
can be seen, as well as its break down at low metallicities of [Fe/H] < —0.7. The Jiang et al. 
(2003) metallicities {magenta squares) are for older clusters in M31 and were derived from 
photometric fits to BC96 models using only three metaUicities, [Fe/H] = 0.0, —0.7, —1.7. 
Their discrete metaUicity steps and any differences between the 1996 and 2003 models are 
the reasons for the large dispersion of these points. 
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For the clusters younger than 1 Gyr, our metaUicites from CN spectrum fits show no cor- 
relation with literature metallicity values. We believe this is because continuum information 
is also necessary for determining metallicities of young clusters, particularly since younger 
clusters likely have higher metallicity with more line blanketing. A hint that the continuum 
information might be helping to determine the metallicity of some young clusters can be seen 
in the plots of our full spectrum fits in Figs. 8 and 9. In Figure 8c, if we ignore the points 
below the 1:1 correlation line for the moment, there does seem to be a metallicity trend with 
the literature for many of the other points. All of the clusters that we ignored below the 
line show possible signs of a near infrared excess in their photometry over the model that 
best fits their optical spectra. Some of these spectra are shown in Figure 3. These clusters 
are in the correct age range, 0.1-2 Gyr, to be affected by TP-AGB stars. It is not clear why 
our derived metallicities of these objects appear too low, since it seems that redder spectra 
would be better fit by models with higher metallicities, but maybe it is more an indication 
that these results cannot be trusted. Fits to the continua shape alone (Figure 8e) result in 
much less correlation to the literature, however, 5 of the 6 outliers below the 1:1 correlation 
line in this case are common to those in the full spectrum plot. Therefore, although it seems 
that correct continuum information might allow metallicity estimates to be made for young 
objects, the uncertainty in proper modeling of the continua shapes makes the full spectrum 
fits unreliable. 

We expect our metallicity estimates that are based on broadband photometry to be less 
well constrained than those based on spectra. Three of our comparison literature sources 
make cluster metallicity estimates based on photometry (denoted in our plots by smaller sized 
symbols): the Gohen & Matthews (1994) VJK colors, the Jiang ct al. (2003) BC96 colors, and 
two of the Santos & Piatti (2004) homogenized literature values. Our photometrically derived 
metallicities (Figure 8d) do not correlate with these literature values, except for possibly a 
very weak relation with the VJK colors {green triangles). As expected, our metallicities 
derived from broadband photometry and the continuum shape show no correlation to any 
values from the literature that are based on spectral lines. 

Santos & Piatti (2004) use empirical relationships between the sums of EWs of Balmer 
and metal lines to the literature age and metallicity values to estimate ages and metallicities 
of clusters in the Magellanic Glouds and in the Galaxy. They find that EW(H(J-|-H7-t-H/?) 
and EW(GaK-l-Gband-l-Mg) are both sensitive to age for clusters younger than 10 Gyr, 
while EW(GaK-l-Gband-l-Mg) is sensitive to [Fe/H] only for clusters older than 10 Gyr. 
Neither of these EW sums correlate with the ages of old clusters {tage > 10 Gyr) or with 
the metallicities of young {tage < 10 Gyr) clusters. Our CN spectrum fits are similar to 
this technique in that they use line strengths, but over the entire spectral range of the data. 
Our fits produce metallicities that agree with literature values only for the older clusters. 
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supporting the Santos & Piatti (2004) result, but extending the minimum age from 10 Gyr 
down to ~ 1 Gyr for vahd metaUicity estimates. 

To summarize our results for metaUicity estimates, we agree very well with the literature 
metaUicities that are based on spectra for clusters older than 1 Gyr. Those based on colors 
have much more scatter, which suggests that our CN spectrum fits are more robust in 
estimating [Fe/H] than methods using colors. For clusters younger than 1 Gyr, metaUicity 
is hard to constrain. Our full spectrum fits hint that correct continua shapes would aid in 
metaUicity estimates of younger clusters, but given the uncertainties in modeling this shape, 
we cannot accurately derive metaUicity estimates for clusters younger than ~ 1 Gyr. 

6.4. Using Model Line Indices 

We compare line indices measured on the BC03 model spectra to those measured on the 
globular cluster spectra in Figure 10 for H/3-[MgFe]' and <Fe>-Mg&. Panels a-d show model 
grids calculated at the spectral resolutions of the data, 6 A, 12 A, 16 A, and 23 A respectively. 
Symbols mark the indices that could be measured for the clusters in each resolution group 
(indices that came out negative due to noise or emission lines were ignored). The 6 and 
23 A resolution groups {squares and crosses) are older clusters in M31, while the rest are 
younger clusters in the Magellanic Clouds. The H/3-[MgFe]' index grids clearly show that 
the model indices become degenerate at tage < 100 Myr for —1.0 < [Fe/H] < +0.5, and at 
tage ^ 1 Gyr for [Fe/H] < —2.0. Most of the young MC clusters fall in the degenerate regions 
of the grids. This degeneracy also occurs for the traditional Lick indices (tabulated with the 
BC03 model package) in this region of parameter space, as shown by the grids in Figure 10c. 
This plot only includes points from clusters that have spectra with resolutions near that of 
Lick indices (11 and 12 A). 

Bruzual & Chariot (2003) investigated how well line indices calculated directly from a 
library of their models containing complex stellar populations with a range of star forma- 
tion histories could match galaxy spectra from the Sloan Digital Sky Survey (SDSS). They 
concluded that those models could simultaneously fit observed strengths of B.f3, H^a, H^a, 
[MgFe]', [MgiFe], [MgsFe], and Dn(4000) in high quaUty galaxy spectra with S/N^ed > 30 
per pixel. Their sample of moderately low redshift SDSS galaxies probes a different region of 
age-metaUicity space than does our sample of globular clusters. Local galaxies are generally 
old with considerably higher metaUicity than these clusters, a region of parameter space 
where the models do well. In contrast, the M31 globular clusters are old and metal poor and 
the MC globular clusters are young and metal poor. Bruzual & Chariot note that their SSP 
models using the STELIB/BaSeL 3.1 library do poorer on line strengths for [Fe/H] < —0.7 
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and our index plots add model degeneracy for tage < 100 Myr, both of which are parts of 
parameter space occupied by globular clusters in this sample. 

Figure lOf shows <Fe> vs. Mg h with model grids included for the minimum and 
maximum spectral resolutions. Some of the clusters depart from the model grid such that 
[Mg/Fe] may be enhanced (Maraston et al. 2003). Since we see signs of a-enhancement and 
the BC03 models use scaled solar abundance ratios, the preferred metal index here is [MgFc]', 
which was shown by Thomas et al. (2003) to be sensitive to mctallicity while insensitive to 
a-cnhanced abundance ratios at the Lick/IDS spectral resolution. However, when [MgFe]' is 
plotted against Balmcr indices, which arc affected by a-cnhanccmcnt because of metal lines 
in the index passbands, many of the clusters fall outside of the model grids with lower Balmer 
indices than the models (Figure lOa-d). One might suspect that the lower H/3 indices are 
caused by the lines being partially filled in by emission, however, many of these clusters are 
old ones from M31 {squares and crosses) and unlikely to have emission. The younger MC 
clusters, on the other hand, could be affected by emission. Nevertheless, since our purpose 
in calculating line indices is for comparison of using the model indices to using the entire 
spectrum, and the affected MC clusters fall in the degenerate region of the model grid with 
unusable indices, we do not correct for emission. 

Further possible causes of the B.f3 index discrepancy could be a-enhancement or the 
presence of blue horizontal branch (BHB) stars in the clusters. The M31 clusters do show 
signs of a-enhanccmcnt in Figure lOf. Line indices are changed by a-enhanccmcnt, but 
in the wrong sense to explain the high H/5 in the BC03 scaled solar abundance ratio model 
grids. H/3 increases with [ct/Fe] (Tantalo & Chiosi 2004a,b; Thomas et al. 2003), due to extra 
absorption in the blue pseudo-continuum and some in the central passband. Therefore, if 
the clusters are a-enhanced their H/3 indices should be higher than the model grid. The 
second possible culprit, BHB stars in old clusters with low metallicities, also affect H/3, but 
in the wrong sense. H/? should increase if a BHB exists in the cluster. We cannot explain 
the H/3 behavior of the clusters relative to the models by either of these causes. 

This same situation is seen by Le Borgne et al. (2004) in their Figs. 12 and 13 where 
they calculate Lick indices directly from the PEGASE-HR model spectra, which use the same 
stellar library as the BC03 models, and compare to globular clusters of intermediate to old 
ages in Andromeda, M31, M33, M81, and M87. They claim that this behavior of H/? is due 
to the clusters being extremely metal-poor, a regime where the interpolated stellar libraries 
suffer from large uncertainties due to lack of sufficient numbers of stars. This same regime 
of [Fe/H]< —0.7 was noted as a less reliable region by Bruzual & Chariot. The metallicites 
of our M31 clusters have literature values of — 1.35<[Fe/H]<0.29, with an average of —0.73, 
so many do fall in the unrehable metaUicity regime of the models. That, coupled with the 
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index degeneracy below ages of 100 Myr, make indices calculated directly from the models 
unusable for most of the globular clusters in our sample. 

Because calculating line indices from the models is troublesome for the age-metallicity 
parameter space occupied by many of the clusters in this sample, we take a different approach 
to investigating the importance of these spectral features. We fit model spectra to the 
clusters only over regions that include and exclude some of the typically used line indices 
and breaks. These features are i:'„(4000), H(5, G4300, H7, H/?, Mg b, Fe5270, and Fe5335. 
The index definitions for the extents of the pseudocontimia passbands of lines, or for the 
flux ratio windows for breaks, define our wavelength ranges over which to fit the models to 
cluster spectra. Figure 11 shows the results of these fits, where a&b use only the regions 
including these features and c&d exclude these regions. It is clear that these features are 
very important for deriving cluster parameters from the models, but also that fitting the 
entire spectral range (refer to Figure 5 and 8) does a better overall job in estimating age and 
metallicity. 



6.5. S/N and Spectral Resolution 

Spectra of z~l galaxies from large surveys will most certainly have lower S/N than the 
spectra of the globular clusters used in the present study. We determine the lowest S/N 
value that can rehably be used to estimate ages by adding gaussian noise of different levels 
to cluster spectra via the noao.artdata.mknoise task in IRAF and fitting CN model spectra 
to them over the wavelength range of 3325-9000 A. We first use actual cluster spectra to 
make the test realistic and avoid any biases that might be introduced by model spectra. The 
clusters in this test were chosen for their high S/N spectra and availability of both age and 
metallicity estimates in the literature. Further tests including effects of spectral resolution 
were then conducted on smoothed model spectra with added noise. 

The resulting ages and metallicities as a function of S/N A^^ for each of four clusters, 
G158, NGC 419, NGC 2134, and NGC 1818, are shown by different symbols in Figure 12. 
The average spectral resolutions of these data are 23A, 13A, 13A, and 13A, respectively. The 
derived ages of the clusters hover around the same values as S/N decreases until ~5 A"^, 
below which the derived ages become unstable. The horizontal fine associated with each set 
of symbols is the average of the derived ages or metaUicites for that cluster from its S/N>5 
spectra. For comparison, the spread of literature ages and metallicities for each cluster are 
shown as vertical lines of the same styles on the right-hand side of the plots. Our derived 
parameters for these four clusters fall within the literature age ranges. If the 0.45 dex [Fe/H] 
scale offset is subtracted from the MC cluster literature values, then our average [Fe/H] lines 
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fall within the hterature values for NGC 419 and G158, which are the clusters older than 
1 Gyr. We do not match the literature [Fe/H] of the younger clusters. We find that the 
spectra need to have S/N > 5 A^^ to achieve stable age estimates from CN spectral fits. 
Examples of fits to NGC 419 with added noise and S/N = 99, 22, and 12 RE-^ (27, 6, 3 A'^) 
are shown in Figure 13. Note that the last one is below the limiting S/N value. 

We have also begun a study of the effects of spectral resolution on derived ages and 
metallicities. To accomplish this we smooth model spectra to AA= 5, 10, 15, 20, 25, and 30 A, 
add noise in 10 realizations at each S/N level up to 60 A^^ for a total of 300-500 simulations 
at each AA, and fit CN model spectra to them. The lowest resolution of AA=30 A was chosen 
to match the lowest typically achieved in galaxy surveys. We select two input models that are 
consistent with age and metallicity combinations occupied by some of the globular clusters 
in this sample and that may also be useful for high redshift galaxies. These combinations 
are: 1 Gyr, [Fe/H]=0.0; and 10 Gyr, [Fe/H]=0.0 (similar to G158 in Figure 12 at 18 Gyr, 
[Fe/H]~ -0.1, AA=23 A). 

The results of these fits are presented in Figures 14-16. Plots in Figure 14 show the 
best fit ages (top rows) and metallicities (bottom rows) as a function of S/N A~^ for spectral 
resolutions of AA = 5, 15, 25 A (to conserve space AA = 10, 20, and 30 A are not shown) 
for the land 10 Gyr, [Fc/H] = inputs. Input parameters arc marked by solid lines, dashed 
lines represent uncertainties of 10% in age and 0.1 dcx in [Fc/H], and dotted lines mark 
age uncertainties of 20%. Small filled circles show individual simulations, open circles arc 
the average derived parameters in S/N bins of width 2, and vertical cyan bars mark the la 
variation within each bin. 

Age-metallicity degeneracy in the derived parameters can clearly be seen in Figure 15 
for the two input models at the same spectral resolutions presented in Figure 14. The 
increased difficulty of distinguishing older ages is also apparent from the much larger spread 
in derived ages for this regime where age changes result in small spectral differences. As 
spectral resolution degrades, derived ages and metallicities for the lower S/N spectra tend 
to be driven to the model grid limits more frequently. 

We choose limiting S/N values for each spectral resolution corresponding to the point 
at which the standard deviation in derived ages for a S/N bin begins to increase beyond 10 
or 20% of the input model age. Figure 16 summarizes the required S/N A~^ as a function of 
AA for the 1 Gyr case on the left and the 10 Gyr case on the right. Solid fines denote derived 
ages within 10% and dashed lines within 20% of the input. Recall that the S/N is calculated 
for these spectra by taking the standard deviation of the fiux in wavelength bins along the 
spectrum, which has the disadvantage that spectral features increase the apparent "noise" 
in the spectrum. Because of this fact, there is a limit to the maximum S/N calculated in this 
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manner. These artificial limits are plotted as dotted lines in Figure 16, and therefore the 
open upward facing triangles are lower limits to the required S/N at those resolutions. The 
10 Gyr object requires higher S/N than the 1 Gyr object at almost all resolutions, implying 
a strong dependence on the object's age. However, because the age is unknown a priori, the 
10 Gyr model must be used to determine the minimum required S/N. 

The expected trend is seen for ages within 10% of the input: higher S/N is required 
for lower spectral resolution. At the 20% level, however, somewhat puzzling behavior is 
seen for the 10 Gyr old object. At AA > 15 A the required S/N actually decreases. This 
does not appear to be an aritfact of the simulations, as we see no correlations between the 
number of points within S/N bins and the standard deviation in derived ages. Interestingly, 
the globular cluster G158 at A A = 23 A in Figure 12 may demonstrate similar behavior. 
Although only one noise realization was performed at each S/N level there, its correct age 
is obtained at much lower S/N than the younger clusters in the plot. Perhaps this trend is 
real. One could envision a scenario in which with decreasing resolution model mismatches 
of narrow lines in the spectrum become less important than broad features typical of old 
spectra, such as the 4000 A break and the Mg absorption trough. In order to get the correct 
age, high S/N is required, but once the allowed uncertainty is increased, the broader features 
can provide enough information in lower S/N spectra. Clearly, this is only an initial guide 
for required S/N and more investigation needs to be done into the cause of this unexpected 
feature. At the high resolution end, our results are consistent with those of Mathis et al. 
(2006) who also fit continuum-removed spectra (via a data compression algorithm) and find 
required S/N~ll-17 A-^ (20-30 per 3 A pixel) at a spectral resolution of 3 A. Extrapolating 
our 5 A point down to 3 A would give S/N~7-22 A"-*^. 

In general, at the spectral resolutions of the globular clusters in Figure 12, AA ~ 13- 
23 A, our simulations indicate that a S/N > 25 A ^ is required for a 20% age uncertainty 
and > 45 A~^ for a 10% uncertainty. This is higher than the S/N~5 A~^ obtained from 
adding noise to actual cluster spectra. This discrepancy is not currently understood and will 
be investigated further in future work. 

7. Discussion 

7.1. Different Spectral Regions 

In §6.4 wc looked at the results of fits to specific regions of the spectrum around typically 
used line indices. Here we test the broader cfi^ects of the regions blueward and redward of the 
4000 A break on derived ages and metallicities. This will be important for fitting galaxies at 
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varying redshifts, since at z~l optical spectra typically reach just beyond the 4000 A break in 
the galaxy's rest frame. For the clusters in our sample that include near ultraviolet spectra, 
we analyze the range of 3200-5650 A, which is the maximum extent common to all clusters in 
this group. Within this range, we fit models to the blue (A < 4000 A), the red (A > 4000 A), 
and the region right around the break (3750-4250 A, as defined in the empirical calibration of 
the break by Gorgas et al. (1999)). The clusters in this study are listed in Table 5 with their 
ages derived from different parts of the spectrum. Correlations of ages from the different 
spectral regions are plotted in Figure 17. 

Fits to different regions give different results. For CN spectrum fits, the red part of the 
spectrum drives the result. The red ages correlate with the ages from the entire spectrum 
and are slightly older than the values in the literature. The blue region gives younger ages 
than both the red region and the literature. Ages from the break region correlate to the red 
ages with about the same amount of scatter as we saw in the red-entire spectrum correlation. 
Since the continuum is normalized out for these fits, the break region is left with only the 
signatures of strong features like the Ca H and K lines and US, similar to the D„(4000)-H(5 
analysis done by Kauffmann et al. (2003a) on SDSS galaxies. The rough agreement of ages 
from this region with those from the entire spectrum emphasizes the importance for the 
4000 A break to be in the spectrum used for estimating the age of an object. The only 
correlation between [Fe/H] is a weak one between the entire spectrum and the red spectrum 
(not shown in Figure 17), not surprising since most metallicity sensitive lines typically used 
are in the red part of the spectrum. 

For full spectrum fits, the trends change and the blue region carries more weight, demon- 
strating the importance of information in the continuum shape of the blue part of the spec- 
trum. This is especially true for young objects with strong Balmer discontinuities. Ages from 
the blue and red regions agree, but with some scatter. However, now the entire spectrum 
agrees equally well with the blue and the red regions. The 4000 A break ages agree much 
better with the blue than with the red ages. The red, blue, and break ages all correlate 
with the literature ages, but the correlation is tightest for the break region, even more so 
than for the full spectrum fits over the entire range. There is no correlation between our 
[Fe/H] derived from the different regions, but the red [Fe/H] has a weak correlation with the 
literature. 

Because of the importance of the continuum information blueward of the 4000 A break, 
we also test hybrid fits that use the full spectrum over A < 4250 A and the CN spectrum over 
A > 3750 A. Wavelengths below 4000 A will not be affected by the TP-AGB and carbon stars 
that dominate the near infrared, so as long as there are no continuum problems with the 
data this type of fit should be vahd. We expected this approach to do better on very young 
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objects. However, the hybrid fits match the full spectrum fits that used the entire wavelength 
range. The continuum shape dominates the fits and spreads out the age correlation with the 
literature for objects younger than ~0.1 Gyr. This is the same problem seen earlier with 
the continua of very young clusters being redder than the models, causing older best-fitting 
ages. 

In summary, the continuum shape contains a lot of age information. In the blue part of 
the spectrum the contimmm is more important than the lines for deriving age, while in the 
red part of the spectrum the lines alone are sufficient to derive the age. The most robust age 
estimates are obtained by using the entire CN spectrum of an object (~ 3200-8200 A), but 
the 4000 A break is an essential feature to be contained within that spectrum. When using 
CN spectra to avoid problems in model and data continua, ages from only the blue region 
are younger than those obtained by the red region or the entire range. Age biases introduced 
by using different regions of the spectrum must be considered when combining objects that 
have different spectral coverage or different redshifts that change the rest frame wavelength 
coverage significantly. Additionally, for objects that cannot be considered a single stellar 
population, the ages determined from different spectral ranges may not be straighforward to 
interpret. Nonetheless, near infrared spectra of high redshift galaxies will be important to 
provide rest frame wavelengths around and redward of the 4000 A break, particularly when 
fitting their CN spectra. 



7.2. TP-AGB and Carbon Stars 

The effects of TP-AGB stars on the near infrared luminosity of objects in which they 
reside have long been known and discussed by groups creating stellar population synthesis 
models (Renzini 1992; Bruzual & Chariot 1993; Girardi et al. 1998; Maraston 1998). Recent 
models by Maraston (2005) incorporate TP-AGB stars through an empirical spectral library 
of C- and 0-type stars by Langon & Mouhcine (2002). The effects of these stars have been 
reanalyzed and cahbrated on MC cluster colors in the new models. They show that over 
80% of the K band light can be contributed by AGB stars between the ages of 0.1 and 2 Gyr 
for solar and lower metallicities. This effect is down to 20% in the V band. 

Maraston fits photometric SEDs of these models to MC cluster SEDs. Two objects 
in common with our sample, NGC 1783 and NGC 419, show discrepancies between the 
Maraston models that best fit the photometric SEDs and BC03 models that best fit the 
spectra. These inconsistencies are the very point being address by the new models, which 
show that younger ages may be obtained for objects that have a near infrared excess due 
to TP-AGB stars if these stars are properly included. However, analysis of the spectra may 



-27- 



reveal a different story. The CN spectra of the two clusters (shown in Figure 18b) are very 
similar to each other. NGC 1783 has D„(4000)=1.295 and NGC 419 has D„(4000)=1.246. We 
derive similar ages for the two clusters by fitting their CN spectra: 1.5 Gyr and [Fe/H]=— 0.8 
for NGC 1783, and 1.5 Gyr and [Fe/H]=-1.2 for NGC 419. The Maraston results for 
NGC 1783 are 0.3 Gyr and [Fe/H] = -0.33, and for NGC 419 are 1 Gyr and [Fe/H]=-1.35. 
BC03 models are overplotted on the CN spectra in Figure 18c, d for all the best- fitting models. 
The older age of 1.5 Gyr is a better fit to the CN spectra of both clusters. 

Our photometric fit age for NGC 1783 is 2.0 Gyr and [Fc/H] = -0.6, and for NGC 419 
is 20 Gyr and [Fc/H] = — 1.4. This much older photometric age for NGC 419 could be due to 
the presence of TP-AGB stars making the continuum redder. The near infrared photometry 
of NGC 419 (see Figure 3) does show an excess over best-fitting models to the spectra. It is 
likely that the BC03 models do not include the enough TP-AGB stars in this case. On the 
other hand, the presence of AGB stars was detected in NGC 1783 by Frogel et al. (1990), but 
apparently their inclusion in the BC03 models is sufficient for this case. We conclude that 
the detailed information contained in spectra is necessary for breaking degeneracies that can 
occur when using only broadband photometry through the near infrared. In a more recent 
paper Maraston et al. (2006) apply the models with enhanced TP-AGB stars to high redshift 
galaxies and find that some degeneracies, particularly between TP-AGB stars and reddening 
due to dust, can be broken when mid-infrared photometry is added from Spitzer. These new 
models look very promising for young objects if infrared data are available, however their 
spectral resolution is not high enough to fit to galaxy spectra. 



7.3. Implications for Galatxies 

The ages expected for early-type galaxies at z~l fall within the range spanned by the 
globular clusters in this sample. The youngest early-type galaxies at these redshifts will likely 
be field galaxies that were the products of mergers. If the bulk of stars in these galaxies 
were formed in mergers at redshifts of z~1.5-2, which are the formation redshifts found for 
field ellipticals (Im et al. 2002; Gebhardt et al. 2003), they would be at least 1.5 Gyr old by 
z~l (for Qm = 0.3, Q\=0.7, Ho=70 km s~^Mpc~^). Our CN spectrum fits estimate ages to 
within 0.16 dex (35%) of the literature values for globular clusters older than 1 Gyr in our 
sample. 

The metallicities of z~l galaxies will likely be higher than many of the clusters in 
this sample, however, a few clusters spanned most of the range expected for galaxies. At 
redshifts of 0.3<z<1.0, galaxies in GOODS-N have metallicities from 0.3 to 2.5 times solar, or 
-0.5 < [Fe/H] < +0A (Kobulnicky & Kewley 2004). We were able to determine the [Fe/H] of 
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globular clusters older than 1 Gyr to within 0.12 dex over the range of —1.6 < [Fe/H] < +0.3. 
The only region we have not tested is the highest metaUicities of [Fe/H] > +0.3. 

Spectra at a resolution of AA = 15 A (R ~ 400 at 6000 A) must have S/N ~ 30- 
45 A""^ in order to derive ages with 20 or 10% uncertainties, respectively, for old galaxies 
and S/N ~ 10-20 A"^ for young galaxies. This is higher than typically achieved for many 
objects in modern galaxy surveys such as MUNICS, DEEP, and TKRS. Our initial checks on 
TKRS spectra from DEIMOS on Keck of GOODS-N galaxies at z~0. 8-1.0 with magnitudes 
of R=23-23.5 find S/N ~ 5 A~^ for individual objects. Therefore, most galaxies will require 
coaddition of a few similar objects to reach sufficient S/N. 

For galaxies at z^l typical optical spectra covering ~4000-9000 A will provide rest 
frame wavelengths of ~2000-4500 A. Our tests have shown that when using CN spectra, 
data to longer wavelengths will provide more robust age constraints. This means that near 
infrared spectra of galaxies will be important data to obtain for accurate age estimates. 

Many early-type galaxies appear to have a-element abundance ratios enhanced above 
solar (Trager et al. 2000a,b; Thomas et al. 2003, 2005). This enhancement is a sign of short 
duration bursts of star formation that may well be typical for early-type galaxies formed 
during mergers, as we expect to find at z~l. We see indications of a-enhancement in 
some of the M31 globular clusters in this sample and we still derive ages within 35% of the 
hterature values when fitting the entire GN spectra. This suggests that results from the entire 
continuum-normalized spectrum, effectively averaging over all the lines, may be less affected 
by a-enhancement than detailed line index methods that rely on only a few affected lines in 
the spectrum. Gassisi et al. (2004) find that broadband colors in the bhie are significantly 
affected by a-enhancement, but removal of the continuum shape must mitigate these effects 
to some extent. 

Model fitting techniques will have to be modified for galaxies. Simple stellar populations 
will no longer be sufficient descriptions of these complex systems. Fits will require models 
with complex star formation histories, multiple stellar populations, and dust. However, this 
work has verified that the reliable regions of model parameter space cover most of that 
expected to be occupied by galaxies. We have also proven techniques by which model fits to 
galaxy spectra can be expected to work. 

8. Summciry and Conclusions 

We have investigated the utility of using the Bruzual & Gharlot (2003) high resolution 
simple stellar population synthesis models for estimating the ages and metaUicities of globular 
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clusters. We chose clusters from M31 and the Magellanic Clouds that cover a wide range 
of age and metallicity values, ~ 0.004<tage(Gyr)<20 and — 1.6<[Fe/H]<+0.3, to explore 
model performance in different regions of parameter space. We tested various techniques of 
fitting the models to data using spectra with and without their continua and photometry. 
Results were compared to age and metallicity estimates in the literature from other methods 
to evaluate our fitting techniques. 

The models are reliable in regions of parameter space where tage >0.1 Gyr and [Fe/H] > 
—0.7. These regions are limited by the incompleteness of the stellar libraries used in the 
models. For objects between the ages of 0.1 and 2.5 Gyr, uncertainties in the amount of near 
infrared light coming from TP-AGB and carbon stars make it hard to match the continuum 
shape of individual globular clusters. Red supergiant stars affect clusters younger than 0.1 
Gyr. Continuum-normalized spectra can be used to get around the TP-AGB star problem, 
as well as any flux calibration issues in the data, with the requirement that high spectral 
resolution models are used. The red supergiant star problem limits accurate age estimates 
to objects older than ~0.1 Gyr. 

Ages from our CN spectrum model fits best match the entire range of age estimates from 
the hterature (Figure 5 and Table 2), with average differences of 0.16 dex (35%) for clusters 
older than 1 Gyr and 0.3 dex (63%) for younger clusters. The methods that depend on the 
continuum shape - the full spectrum, continuum, and photometry fits - overestimate the 
ages of many clusters below 1 Gyr. Wc believe this to be due to a near infrared excess from 
TP-AGB stars making the cluster spectra redder than the models and forcing older aged 
models to provide the best fits for cluster ages of 0.1-1 Gyr. Red supergiant stars similarly 
affect clusters younger than 0.1 Gyr. Although photometric fits seem to do well for clusters 
older than 1 Gyr (0.1 dex or 23% age error relative to the hterature), the CN spectrum fits 
are better for simultaneously providing both age and metaUicity for these older objects. Fits 
to photometry alone cannot reproduce the correct metallicities. 

Comparing the fitting methods (Figure 4), we find a tighter correlation between ages 
derived from the full spectrum and continuum fits than from the full spectrum and CN 
spectrum fits for clusters younger than 1 Gyr, suggesting that these derived ages are more 
strongly driven by the continuum shape than by the spectral lines. Clusters older than 1 Gyr 
show a tighter correlation between ages derived from the CN spectrum and photometry fits. 
Because the derived ages of young clusters are so strongly influenced by the continua and we 
have seen signs of problems in the models matching the cluster continua, wc conclude that the 
best method for deriving accurate ages for globular clusters, especially when simultaneously 
determining metallicity, is fitting models to their continuum-normalized spectra. 

On metallicities of clusters older than 1 Gyr, we agree to 0.12 dex with the literature 
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values that are based on spectra. Those based on colors have much more scatter, which 
suggests that CN spectrum fits are more robust in estimating [Fe/H] than methods using 
colors. We cannot accurately derive metallicity estimates for clusters younger than ~ 1 Gyr. 

Blueward of the 4000 A break the continuum shape is more important than the lines 
for deriving age. Redward of the break the lines alone in a continuum-normalized spectrum 
are sufficient to derive the age. Many of the globular clusters tested here have spectra out 
to rest wavelengths of 1 ^m. If fitting models to continuum-normalized spectra of z~l 
galaxies, this means that near infrared spectra through K band will be important to provide 
similar spectral coverage and constraints on age and metallicity. The 4000 A break is a 
very important feature to be contained in the spectrum. It will be contained in typical z~l 
spectra, but higher redshifts may require near infrared spectra to correctly constrain ages 
and metaUicites. We should note, however, that although a strong age indicator for globular 
clusters and early type galaxies, the 4000 A break can be significantly affected by dust in 
star forming galaxies. For example, MacArthur (2005) find dust imparted age errors as large 
as 9.5 Gyr for a 0.5 Gyr BC03 model SSP when ages are derived from the D„(4000)-Fe4668 
plane. 

The spectral S/N required to derive ages by fitting the entire CN spectrum is similar 
to the 30-40 A^^ typically required when using line indices on old objects. For a 10 Gyr 
old object at AA = 5-25 A our required S/N is 15-30 A~^ for a 20% age uncertainty and 
30-55 A~^ for a 10% uncertainty. A 1 Gyr year old object requires S/N ~ 3-15 A^^ for a 
20% age uncertainty and 10-25 A^-*^ for a 10% uncertainty. At AA = 30 A both ages require 
S/N > 60 A. Fairly low resolution, R~240-1200, spectra can be used, at the cost of increasing 
required S/N with decreasing resolution. 

We have successfully tested most of the age-metallicity space expected to be occupied by 
galaxies at z~l, except for the highest metallicities of [Fe/H]>-|-0.3. Fitting the continuum- 
normalized spectra of these galaxies should work well for spectra without accurate flux 
calibration, typical when faint galaxies are observed on multiple masks over long periods of 
time during surveys. The required S/N and spectral resolution for determining ages and 
metallicities matches that typically achieved in recent galaxy surveys, as long as a few like 
objects arc coadded to boost the S/N. The techniques presented here will be apphed to 
galaxies in future work. 

We wish to thank the anonymous referee for very useful comments and suggestions. 
MJW was partially funded for this work by NASA GSRP grant number NGT 5-50301. 
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Fig. 1. — The effect of S/N on Ax^ contours of age and metallicity for CN spectrum fits to 
tlie SMC globular cluster NGC 419. The S/N per resolution element from top left to lower 
right is 99, 22, 12, 8, 6, and 4. The best-fitting age and [Fe/H] is given in the plot titles and 
marked by the black crosses. The yellow circles mark the age and [Fe/H] from the highest 
S/N fit. The solid, dashed, and dotted lines are 1, 2, and 3a error contours, respectively. 
Color bars give the Ax^ values for 68.3% (Ax^=2.3), 90%, 95.4%, 99%, 99.73%, and 99.99% 
(Ax^=18.4) confidence levels for two degrees of freedom. (The highest level is converted to 
white for plot clarity.) Black dots mark the locations of models in the age-metallicity grid. 
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Fig. 2. — Example model fits to two clusters of very different ages, G177 in M31 and NGC 
1711 in the LMC a & b are the full spectrum fits, c & d are the continuum-normalized 
spectrum fits, and e & f are the continuum fits. Solid black lines are the cluster spectra, 
dashed cyan lines are model spectra corresponding to the best spectroscopic fit, dotted 
green lines are the model spectra corresponding to the best photometric fit, filled circles 
are the cluster photometry, triangles are photometry from the best photometric fit models, 
and horizontal bars mark the filter bandwidths. Photometry for G177 includes U,B,V,J,K 
bands and NGC 1711 includes U,B,V,J,H,K bands. Best-fitting ages and metallicities for 
each method are given in the plot labels. In the continuum fits, black lines are the cluster 
continua and green lines are the best-fitting model continua. 
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Fig. 3. — Spectra with possible signs of TP-AGB or carbon stars because their near infrared 
photometry points show an excess over the best fitting model to the spectrum. The legend is 
the same as in Figure 2. The ages of these clusters all fall within the range where TP-AGB 
and carbon stars are important, 0.1-2 Gyr. Our CN spectrum fits give ages of 0.1 Gyr for 
NGC 1850, 1.5 Gyr for NGC 419, 0.2 Gyr for NGC 2134, and 0.1 Gyr for NGC 2136. 
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Fig. 4. — Comparison of ages and metallicities derived by different types of model fits for 
LMC clusters [crosses), SMC clusters (circles), and M31 clusters (asterisks) with linear fits 
(solid red lines) and lines of 1:1 correlation (dotted lines). Clusters with average literature 
ages > 1 Gyr are marked by yellow circles and < 1 Gyr by cyan diamonds. The error bars 
are la. If no error bars are seen, they are smaller than the symbols. Age comparisons are 
between a) the full and CN spectrum fits, b) the photometry and CN spectrum fits, c) the 
continuum and full spectrum fits, and d) the continuum and photometry fits. 
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Fig. 5. — Cluster ages derived from model fits compared to ages from the literature, a) 
is the CN spectrum fits, b) is the full spectrum fits, c) is the photometry fits, and d) 
is the continuum fits. Literature ages are derived from main sequence turnoffs by Hodge 
(1983) {black crosses), UBV photometry by Santos et al. (1995) {blue asterisks), CMDs by 
Girardi et al. (1995, 1998) and Elson & Fall (1988) {green diamonds), photmetric fits to 
BC96 models by Jiang et al. (2003) {magenta squares), line index ratios by Leonardi & Rose 
(2003) {red circles), and literature ages on a homogeneous scale by Santos & Piatti (2004) 
{black inverted triangles). The thick horizontal cyan lines connect points that are the same 
cluster with multiple literature ages. Thin horizontal lines are quoted error bars on literature 
ages, thin vertical lines are la error bars on our derived ages. If no error bars are seen, they 
are smaller than the symbols. 
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Fig. 6. — Absolute values of fractional age errors relative to literature values for CN spectrum, 
full spectrum, photometry, and continuum fits. The dotted lines mark the 100% error point. 
Symbols are the same as in Figure 5 with literature ages from main sequence turnoffs by 
Hodge (1983) {black crosses), UBV photometry by Santos et al. (1995) {blue asterisks) and 
Girardi et al. (1995, 1998) {green diamonds), photmetric fits to BC96 models by Jiang et 
al. (2003) {magenta squares), line index ratios by Leonardi & Rose (2003) {red circles), and 
literature ages on a homogeneous scale by Santos & Piatti (2004) {black inverted triangles). 
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Fig. 7. — Spectra of extreme outliers in the age correlation plot for CN spectrum fits in 
Figure 5. Here, a&b) are spectra of NGC 2092 and c&d) are spectra of NGC 2096. a&c) are 
the full spectrum fits to show the original spectral shapes, and b&d) are the CN spectrum 
fits. After masking emission lines, NGC 2092 does not have enough distinguishing spectral 
features to obtain a good fit. NGC 2096 appears older than its literature age of 49 Myr from 
its 4000 A break strength of D„(4000)=1.28. These clusters are both known to contain red 
supergiants and NGC 2092 is embedded in a nebula. 
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Fig. 8. — Cluster [Fe/H] derived from our different types of model fits compared to the 
literature, a&b) are from the CN spectrum fits for old and young clusters separately, c) is 
from the full spectrum fits, d) is from the photometry fits, and e) is from the continuum 
fits. Literature values are from comparison values in Leonardi & Rose (2003) (which come 
from Olszewski et al. (1991), Cohen (1982), Seggewiss & Richtler (1989), Sagar & Pandey 
(1989), and Piatti et al. (2002)) in the MC {black crosses), line indices by Huchra et al. 
(1991) in M31 {blue pluses), line index ratios by Leonardi & Rose (2003) in the MC {red 
circles), VJK colors by Cohen & Matthews (1994) in M31 {green triangles), photmetric fits 
to BC96 models by Jiang et al. (2003) in M31 {magenta squares), and literature [Fe/H] on 
a homogeneous scale by Santos & Piatti (2004) in the MC {black inverted triangles). Large 
symbols denote literature values based on spectra and small symbols are values based on 
colors. The thick horizontal cyan lines connect points that are the same cluster with multiple 
literature metallicities. Thin horizontal lines are error bars on literature [Fe/H], thin vertical 
lines are la error bars on our derived values. The best agreement with literature metallicity 
occurs for spectra-based values of clusters older than 1 Gyr from our CN spectrum fits. 
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Fig. 9. — [Fe/H] errors relative to literature values (fit - literature) for a) CN spectrum, b) 
full spectrum, c) photometry, and d) continuum fits, plotted against average literature ages. 
Symbols are the same as in Figure 8 with [Fe/H] from comparison values in Leonardi & Rose 
(2003) in the MC {black crosses), line indices by Huchra et al. (1991) in M31 {blue pluses), 
line index ratios by Leonardi & Rose (2003) in the MC {red circles), VJK colors by Cohen 
& Matthews (1994) in M31 {green triangles), photmetric fits to BC96 models by Jiang et al. 
(2003) in M31 {magenta squares), and literature [Fe/H] on a homogeneous scale by Santos 
& Piatti (2004) in the MC {black inverted triangles). Large symbols denote literature [Fe/H] 
derived from spectra and small symbols are from colors. The dependence of [Fe/H] errors 
on cluster age can more easily be seen here for all types of fits. 



-48- 









1 £ ciLElers " 






















■ 










-1.0 — -^-^jo 




■ QO^^ 













2 

IMgFe]' 






33 ciLElers " 














« 1 , 1 . 






Fig. 10. — Line indices calculated from BC03 models, a-d) H/9-[MgFe]' model indices calcu- 
lated at different spectral resolutions. Grid age lines (solid) from bottom to top are 20, 10, 5, 
2, 1, 0.1, 0.01, and 0.001 Gyr, and [Fe/H] lines (dashed) from right to left are +0.5, 0.0, -0.5, 
-1.0, and -2.0 (these are labeled in panel a). Note that the model indices become degenerate 
at tage < 100 Myr for -1.0 < [Fe/H] < +0.5, and at tage < 1 Gyr for [Fe/H] < -2.0. Cluster 
indices are plotted as symbols, differientiated by spectral resolution: 6 A (squares, M31), 
11 A (triangles, MC), 12 A (circles, MC), 16 A (diamonds, MC), and 23 A (crosses, M31). 
e) contains grids of the traditional Lick indices provided with the BC03 model package. 
They suffer from similar degeneracy at low ages and metallicities. f) shows an indication of 
[Mg/Fe] enhancement in some of the M31 clusters. Model grids for indices at the minimum 
and maximum spectral resolutions are included. 
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Fig. 11. — Ages and metallicities from CN spectrum fits to only the regions around typically 
used line indices and spectral features (a&b) and to the regions excluding these features 
(c&d). The features used were D„(4000), Rd, G4300, H7, H/3, Mg b, Fe5270, and Fe5335. 
Refer to the captions on Figure 5 and 8 for symbol description. 
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Fig. 12. — Dependence of ages and metallicities derived from CN spectrum fits on average 
S/N A^-*^. S/N=5 is marked by vertical lines. The clusters shown are G158 {squares, solid 
red line), NGC 419 {pluses, dashed green line), NGC 2134 {triangles, dotted blue line), and 
NGC 1818 {circles, dash-dot magenta line). The horizontal lines mark the average age or 
metallicity of each cluster for the S/N>5 points. The spectral resolution of these clusters 
is 23A for G158 and 13 A for the others. The vertical bars on the right show the spread of 
literature values for each cluster, in the same line types and colors. If the 0.45 dex [Fe/H] 
scale offset is subtracted from the MC cluster (NGC 419, 2134, 1818) literature values, then 
our derived [Fe/H] for NGC 419 falls within the literature range. Our [Fe/H] values for the 
two younger clusters do not. 
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Fig. 13. — NGC 419 best-fitting spectra with artificially added noise resulting in S/N~99, 
S/N~22, and S/N~12. Between the last two S/N levels is where the model fits begin to give 
unstable ages for this cluster. 
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Fig. 14.— Noise simulations for 1 and 10 Gyr, [Fe/H] = 0.0 at AA = 5, 15, 25 A. The 
top rows are derived age and bottom rows are derived [Fe/H]. Input age and metallicity are 
marked by solid lines, age uncertainties of 10% by dashed lines, age uncertainties of 20% by 
dotted lines, and [Fe/H] uncertainties of 0.1 dex by dashed lines. Small filled circles represent 
individual simulations, open circles mark the average derived parameter in S/N bins of width 
2, and vertical cyan bars mark the la variation in derived values within each bin. 



- 53 - 



1 2 30 ^0 50 60 



10.0 




1.0 
0.5 



0.0 



-0.5 ^ 



-1.0 



dinJaUcTE 

i lll |l II l| III l| II II |l III |l I N 



390 dnUjUcTE 

fr... I I ... .1 .... I. ... I. ..7 



' ^pl: 10(^, F«*fl-O0. il.-1SA 
.1 ... .1 .... I. ... I I ... . 



300 ShlLtaUOTG 

. .1 ... .1 .... I. ... I I ... . 



470 stniJdorE 
rpil: 10 (^r, [J^ifl-D-aO il-SA 



I n il I III I II I III I II II 



470 stnUdlom 
. .1 .... I. ... I I ... .1 .. .r 



lOjQ 



IjO 
0.5 



OJO 



- -0.5 



-1.0 



10 20 30 * 50 60 



1 20 30 ^10 50 00 



Fig. 14. — cont. 
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Fig. 15. — Age-metallicity plots of noise simulations with AA = 5, 15, 25 A. The top row 
is 1 Gyr and bottom row is 10 Gyr. Symbol size denotes the S/N A~^: largest black circles 
have S/N > 20, next smaller cyan circles have 10 < S/N < 20, next smaller red circles have 
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Fig. 16. — S/N required to derive ages within 10% (solid lines) and 20% (dashed lines) 
of a 1 Gyr input model on the left and a 10 Gyr model on the right (both with [Fe/H] = 0) 
for a range of spectral resolutions. The dotted lines mark maximum S/N values at each 
resolution (see text). The upward facing open triangles are lower limits to the required S/N. 
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Fig. 17.— Age results from fitting the blue (A < 4000 A), red (A > 4000 A), and 4000 A break 
(3750-4250 A) regions of the spectra, a-e) are CN spectrum fits, f-j) are full spectrum fits, 
k-n) are CN spectrum fit ages compared to the literature, and o-r) are full spectrum fit ages 
compared to the literature, h&r) are from using the entire spectral range. 
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Fig. 17. — cont. 
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Fig. 18. — NGC 1783 and NGC 419 were two clusters in common with Maraston (2005) who 
fit enhanced TP-AGB star models to the photometric SEDs. a) shows the two cluster spectra 
from Santos et al. (2002), b) shows their CN spectra, and c&d) include BC03 model CN 
spectra. Our best fitting model to NGC 1783 was 1.5 Gyr and [Fe/H]=-0.8 and to NGC 419 
was 1.5 Gyr and [Fe/H] =— 1.2 {blue dotted lines). Maraston's best fit to the photometry of 
NGC 1783 was 0.3 Gyr and [Fe/H]=— 0.33 (the nearest model in our grid was 0.4 Gyr and 
[Fe/H] = -0.3) and of NGC 419 was 1 Gyr and [Fe/H]=-1.35 {red dot-dashed lines). 
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Table 1. Cluster Wavelength Coverage and Resolution 





UV 


UV 


UV 


Optical 


Optical 


Optical 


IR 


IR 


IR 




Object 


^start 


^end 


Resolution 


^start 


^end 


Resolution 


^start 


^end 


Resolution 


Refs 


LMC 


H88-267 








3500 
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3500 


5870 


16 








1 


HD32228 
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Tablc 1 — Continued 



UV UV UV Optical Optical Optical IR IR IR 

Object Xstart ^end Resolution Xstart ^end Resolution Xstart X^nd Resolution Refs 



NGC1984 








3500 
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NGC1994 
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NGC2002 
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5656 
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NGC2009 
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NGC2011 
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14 


1,3 


NGC2100 
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14 
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16 








1 
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4704 


12 








1 
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16 








1 
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11 








5 


NGC299 
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16 


5600 


10000 


14 
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5656 
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11 


5658 


10116 


14 


2,1,3,5 
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5654 


15 


3700 


8200 


11 








2,5 



M31 
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Boll 18 



3200 6100 6 

3200 6100 6 



7 
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Table 1 — Continued 



UV UV UV Optical Optical Optical IR IR IR 

Object Xstart ^end Resolution Xstart X^nd Resolution Xstart Xf.„d Resolution Refs 



B01119 


3200 


6100 


6 


7 


B01124 


3200 


6100 


6 


7 


Gl 


3600 


9900 


23 


8 


G78 


3600 


9900 


23 


8 


G158 


3600 


9900 


23 


7,8 


G165 


3200 


6100 


6 


7 


G168 


3200 


6100 


6 


7 


G169 


3200 


6100 


6 


7 


G170 


3600 


9900 


23 


8 


G174 


3200 


6100 


6 


7 


G175 


3200 


6100 


6 


7 


G177 


3600 


9900 


23 


7,8 


G184 


3200 


6100 


6 


7 


G185 


3200 


6100 


6 


7 


G190 


3200 


6100 


6 


7 


G222 


3600 


7934 


8 


8 



References. — (1) Santos et al. (1995), (2) Bica et al. (1994), (3) Bica et al. (1990), (4) Bica ct al. (1987a), (5) Bica et 
al. (1986a), (6) unpublished, Santos et al. (2002), (7) Jablonka et al. (1998), (8) Jablonka et al. (1992) 



Table 2. Offsets and Dispersions of Derived Cluster Ages 





(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 






All 






UBV 


Line Index 


Homog. 


T>/^n/^ Til J. 

BC9d Fnot. 


Old 


Young 




All 


Good 


MSTO 


CMDs* 


Colors 


Ratios 


Scale 


Fits 


> 1 Gyr 


< 1 Gyr 


U Ct Lr Ivjilo 




MP M31 


MC 


MC 


MC 


MC 


MC 


M31 


MC M31 


MC 


R.cfoi'onccs 


1.2.3.1.5.6 


1.2. .3. 1.5.6 


1 


3.5 


2 




6 


4 


1.3.4,5.6 


1.2.3.5.6 


No. objects 






29 


11 


51 


11 


18 


8 


1 


(il 


^nf^ci^nim Pit* 
'I'l Ulli -T llio 






















A VP- Off^pt rOvr^ 

jT. V g vyilDCt I vJrj'l J 


0.84 


0.19 


0.12 


0.56 


0.03 


0.26 


0.59 


4.36 


1.63 


0.09 


A VP- OfFciPt Mpv^ 


0.34 


0.29 


0.36 


0.21 


0.29 


0.22 


0.31 


0.25 


0.16 


0.30 


RMS (Gyr) 


3.18 


1.00 


0.82 


1.06 


0.12 


0.79 


1.82 


8.34 


3.78 


0.35 


RMS Cdex'l 


0.57 


0.36 


0.45 


0.31 


0.37 


0.27 


0.40 


0.46 


0.17 


0.37 
























Ave- Offspt rOvr") 


0.93 


0.54 


0.71 


2.54 


0.42 


2.24 


0.70 


3.24 


2.76 


0.39 




0.70 


0.69 


0.64 


0.36 


0.78 


0.33 


0.52 


0.36 


0.70 


0.69 


RMS (Gyr) 


4.42 


3.20 


3.68 


6.08 


1.56 


5.78 


5.43 


9.61 


11.55 


1.47 


RMS (dex) 


0.93 


0.90 


0.71 


0.51 


0.98 


0.44 


0.64 


0.42 


0.79 


0.91 


Photometry Fits 






















Avg Offset (Gyr) 


0.99 


0.48 


0.73 


0.67 


0.34 


0.37 


0.74 


5.24 


2.06 


0.38 


Avg Offset (dex) 


0.004 


0.59 


0.71 


0.42 


0.63 


0.38 


0.41 


0.32 


0.10 


0.62 


RMS (Gyr) 


3.29 


1.19 


1.55 


0.85 


1.04 


0.61 


1.42 


9.47 


2.67 


1.02 


RMS (dex) 


0.52 


0.76 


0.90 


0.52 


0.80 


0.48 


0.49 


0.49 


0.11 


0.78 


Continuum Fits 






















Avg Offset (Gyr) 


0.87 


1.17 


2.05 


4.51 


0.46 


4.21 


2.90 


-2.49 


6.06 


0.85 


Avg Offset (dex) 


0.06 


0.80 


0.84 


0.49 


0.84 


0.46 


0.73 


0.55 


0.80 


0.80 


RMS (Gyr) 


5.44 


4.30 


5.51 


8.66 


1.57 


8.29 


7.71 


10.88 


12.72 


3.02 


RMS (dex) 


0.74 


1.03 


1.05 


0.68 


1.04 


0.58 


0.93 


0.65 


0.73 


1.04 



Note. — Column 1 includes all clusters with literature ages and column 2 has excluded the two outliers in the upper left plot of Figure 5 and the 
objects that only have BC96 photometric fits from the literature. All columns to the right, except for column 8, also exclude these points. 

"■These age estimates were collected from Girardi et al. (1995, 1998) and Elson & Fall (1988) by Leonardi & Rose (2003) for comparison to their 
estimates from line index ratios. 



References. — (1) Hodge (1983), (2) Santos et al. (1995), (3) Girardi et al. (1995, 1998), (4) Jiang et al. (2003), (5) Leonardi & Rose (2003), (6) 
Santos & Piatti (2004) 



Table 3. Offsets and Dispersions of Derived Cluster [Fe/H] 





(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 






Line 


V.JK 


BC96 


LIR 


Line Index 


Hornog. 


Old 


Yonng 




All 


lu(li(.-(*s 


Colors 


( 'olors 


Li! ( 'oiiii)ai(*'^ 


Hatios 


Sc-ale 


1 (tVI- 


< 1 ^^y^ 


l^OCcttlOHS 


iVi^, IVIOJ. 


M31 


M31 


M31 


MC 


MC 


MC 


A/TP A/T'il 


MC 


References 


1 O Q /I A 


1 
1 


n 
A 


q 


A 
4t 


u 


a 
D 


1 O Q /I K 
i,2,3,4,b,b 




No. objects 




a 
D 


Q 


Q 
C5 


1 1 

J.1 


c 




1 8 


1 9 


91 


CN Spectrum Fits 




















Avg Offset (Gyr) 


-0.23 


0.11 


0.25 


-0.14 


-0.68 


-0.69 


-0.37 


-0.10 


-0.29 


RMS (Gyr) 


0.85 


0.12 


0.52 


0.71 


1.06 


0.81 


0.91 


0.41 


1.02 


Full Spectrum Fits 




















Avg Offset (Gyr) 


0.23 


0.31 


0.88 


0.40 


-0.27 


-0.05 


0.16 


0.48 


0.10 


RMS (Gyr) 


0.84 


0.50 


1.03 


0.59 


0.89 


0.53 


0.71 


0.78 


0.88 


Photometry Fits 




















Avg Offset (Gyr) 


0.28 


0.62 


0.77 


0.41 


-0.18 


-0.11 


0.15 


0.51 


0.14 


RMS (Gyr) 


0.87 


0.74 


0.84 


0.67 


1.10 


0.89 


0.72 


0.68 


0.96 


Continuum Fits 




















Avg Offset (Gyr) 


0.36 


0.81 


1.34 


0.85 


-0.29 


-0.35 


0.12 


0.85 


0.08 


RMS (Gyr) 


1.11 


1.03 


1.49 


0.92 


1.02 


0.92 


0.98 


1.09 


1.13 



Note. — Column 1 includes all clusters with literature [Fe/H]. 

'^The Leonardi & Rose (2003) comparison literature values come from Olszewski et al. (1991), Cohen (1982), Seggewiss & Richtler (1989), 
Sagar & Pandey (1989), and Piatti et al. (2002). 



References. — (1) Huchra et al. (1991), (2) Cohen & Matthews (1994), (3) Jiang et al. (2003), (4) Leonardi & Rose (2003) literature 
comparison values, (5) Leonardi & Rose (2003), (6) Santos & Piatti (2004) 



Table 4. Cluster Ages and Metallicities 





CN Spec 


CN Spec 


Full Spec 


VBY^ MSTO'' LIR<^ HLit<^ BC96= Ll' 


VJKS 


LIR lit^ 


LIR'^ 


HLit^ 


Object 


Age 


[Fe/H] 


[Fe/H] 


Age Age Age Age Age [Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 



LMC 


H88-267 


001+" ""^ 

U.UU1_Q QQQ 


-2 20+" S 


-1 80+"°5 


HD32228 


002+° °°^ 


50+° °? 

0.05 


-2 20+° °S 


HD269828 




-0 7o+°?^; 

'"—0.10 


-1 10+°-°? 

.'-•.'-" — 0.05 


IC2105 


U.UUD_Q „Qj 


-2 2o+!;-"^ 

^■^"—0.00 


-2 00+°-?? 

^•""—0.10 


IC2111 


n nnc+O-OOl 


-1 60+°-°? 

-L.uu_Q 05 


-1 90+°/'? 

^•-^"-0.05 


IC2116 




50+° °2 


-1 20+S°? 

^•^"-0.05 


KMHK1074 


004+° °°l 


-0 10+°°? 

"•^^—0.05 


10+005 


NGC1466 


2o.oooi°:°r,S 


-2 00+°-°-^ 

^■"" — 0.05 


-1 30+°-°? 


NGC1711 


0.080«:«1« 


-1 40+of » 

-^-■^"-0.05 


-0 20+!;-"? 

"•^"—0.05 

30+25!? 


NGC1714 


20.000l»:««« 


"•^"—0.05 


NGC1731 





-0.000 


-2 20+° °^ 
^•■^'-'—0.00 


-1 10+°-°? 

^•^"-0.05 


NGC1735 





"^"-0.010 


-2 20+° °? 

^•■^"—0.00 


-1 60+°°? 

^•^'-'—0.05 


NGC1736 







-1 90+°?° 

-^■^"—0.10 


-1 70+°-°? 


NGC1743 





Q1Q + 0.005 

-,(,(.+0.050 
1*^^-0.010 


-2 20+°-»° 

^■^"—0.00 


-0 70+;;-?? 

"• ' "—0.10 


NGC1755 





-1 70+?;°5 

^- '"-0.05 


-1 30+°?? 

^•""-0.10 


NGC1761 





^,-,(,+0.005 
"1^-0.001 


-2.20+°°5 
— u.uu 


0.50+°°5 


NGC1767 





QOg+0.001 
L'LID-O.OOl 


10+°°^ 

"■J^"_0.05 


20+° °^ 
U-^U_o 05 


NGC1772 





(,20+0.010 
'^^^-0.005 


n 50+0 00 

'^•'^'J-0.05 


-0 60+° "'' 

'^•"'^-0.05 


NGC1774 





-,(,(,+0.050 
I'^'^-O.OIO 


-1 90+"-"-''' 
i--'"-o.io 


1 sn+0.05 

-l-«"-0.05 


NGC1782 





(,-(,+0.010 
'-'^'-'-0.010 


-0.401°:- 


-0 30+0 05 
'^••^'^-0.05 


NGC1783 


1 


^-(,(,+0.250 
'^'J^-0.160 


-o.8oin^ 


00+0 05 
U.UU_o 05 


NGC1805 





nns+0 001 


co+o oo 

U-OU_o 05 


Ko+o oo 

U-OU_o.05 


NGC1818 





(,-(,+0.010 
^^^-0.010 


30+°-°-°^ 

l-''"-0.05 


-0 50+0 05 
'^•"'^-0.05 

-1 50+0-05 
l-"U-0.05 


NGC1831 





onn+o-ioo 

*UU_ 0.100 


1 cn + 0.05 


NGC1847 





(,2Q+0.010 
'-'^'-'-0.005 


_-, 70+0.05 

1- "-'-0.05 


10+005 
^•1^-0.05 


NGC1850 





-,(,^,+0.050 
1""-0.010 


-1 60+°°^ 
l.DU_o 05 


-1 60+°-°^ 
l.DU_o 05 


NGC1850A 





002+0.001 
L'LI'^-O.OOl 


c;o+0 00 
U-OU_o JO 


9 ofl+O.OS 
^•'2U_o 00 


NGC1854 





0(50+0.010 
i°i-o.oio 


_1 40+0.05 

l-^"-0.05 


10+"-°'' 


NGC1856 





op,f.+0.100 
0.100 


-I cn+0.05 


-0.901°:- 


NGC1866 





200+0.100 
^^^-0.050 


-2 20+°-"-''' 
^■^'^-0.00 


-0.801°:- 


NGC1868 


3 


000+0-500 
"'-'"-0.600 


_-, -0+0.05 


-0 10+°°5 
U-lU_o 05 


NGC1872 


1 


000+0 100 

L'"L)_0.100 


9 on+o.os 


00+° °^ 

U.UU_o 06 



0.005 
0.004 
0.004 



0.004 

0.020 

0.004 
0.032 

0.055 

0.007 

0.009 
0.037 
0.036 

0.009 
0.014 

0.021 
0.031 
0.004 
0.034 



0.136 



1.200 

017 
190 
016 
040 

030 
120 
080 
330 



.548 

020 
501 



302 
933 



.100 
068 



300 
014 

1.320 

1.031 

034 
151 
150 
850 



1.75 

90 
01 



1.54 



1.65 



20 
50 



1.32 



Table 4 — Continued 





CN Spec 


CN Spec 


Pull Spec 


UBV^ 


MSTO'' 


LIR« 


HLit'^ Bcge" 


Lif 


VJKfs 


LIR lit'' 




HLit^ 


Object 


Age 


[Fe/H] 


[Fe/H] 


Age 


Age 


Age 


Age Age 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 



NGC1895 
NGC1903 
NGC1921 
NGC1935 
NGC1936 
NGC1949 
NGC1951 
NGC1967 
NGC1970 
NGC1978 
NGC1983 
NGC1984 
NGC1994 
NGC2002 
NGC2003 
NGC2004 
NGC2006 
NGC2009 
NGC2011 
NGC2021 
NGC2025 
NGC2032 
NGC2033 
NGC2035 
NGC2041 
NGC2058 
NGC2070 
NGC2091 
NGC2092 
NGC2096 
NGC2098 
NGC2100 
NGC2102 



0.040«:»1° 
0.080l0°l» 
200+°- 

U.ZUU_Q QgQ 

002+° °°^ 
100+0060 

U.1UU_Q QJQ 

n nnc+0.001 

U.UZU_Q QQg 

010+0 005 

U.U1U_Q QQ]^ 

002+0 001 
4 000+0-500 

002+0-001 

0.00410.0^1 

0.00810-021 

0.008lO;Oo°ol 
006+0-001 

010+0-005 

U.UiU_Q QQ]^ 

006+0-001 

U.UUD_Q QQ]^ 

010+0-005 

U.U1U_Q QQ]^ 

006+0-002 

U.UUD_(j QQj^ 

020+0-010 

U.UZU_Q QQg 

100+0 050 

U.lUU_g Qj^Q 

040+0 010 
002+0-001 
002+0-001 
100+0-050 

U.1UU_Q Q^Q 

200+0-100 

U.ZUU_Q 

060+0-010 

U.UDU_Q QJ^Q 

0.020tO;OJ2 
10.000i?:500 

6 000+0-500 

020+0-010 

006+0-001 

004+0-001 



-0 70+0-05 
-2 00+0-05 
-2 20+0-05 

-1+0.40 
-0.10 
-1+0.10 
-0.20 
-1+0.05 
-0.00 
-1+0.10 
-0.10 

-0 30+0-05 

-1+0.05 
-0.00 
-1+0.05 
-0.05 
30+0 05 

0.50^0-00 



-2.103 
-2.00j 
-2.20j 
-1.80"* 



-2.20: 
-1.30"* 



o.3o: 

o.4o: 

o.io; 

o.io; 

o.oo; 

o.3o: 

o.3o: 

-o.8o: 

-i.9o: 

-2.20: 

o.3o: 
-i.9o: 

-1.70 
-1.50 
0.50 
0.10 
-1.70 
-2.20 
0.00 
0.20 
-0.10 



).05 
).()5 
).05 
).05 
0.05 
1.05 
-0.05 
-1+0.05 
-0.05 
-1+0.05 
-0.05 
fO.05 
-0.05 
-1+0.05 
-0.10 
-1+0.10 
-0.05 
fO.lO 
-0.05 
fO.05 
■'-O.OO 
fO.05 
-0.05 
-1+0.40 
-0.30 

-1 70+010 

l--'"-0.05 
+0.10 
0.30 
+0.00 
0.05 
0.05 
0.05 
+0.05 
0.10 
+0.05 
0.00 
+0.05 
0.05 
+0.05 
0.05 
+0.06 
0.05 



-0 90+0-05 
'^•^"^-0.05 

_-, 7Q+0.05 
l-'^-0.05 

-2 10+0-10 

F0.05 



-1.50 



-0.05 
+0.05 



0-20_o.io 

_-, 7Q+0.05 

-0 30+0-05 

'^••^'^-0.05 

-0.8010;- 

hO.05 
-0.05 
0.05 
0.05 
-1 50+0-05 

i-O'J-o.os 

-,+0.05 
-0.05 

_9 r,ri+0.05 

^■^"-0.00 

20+0-05 
'^■^"-0.05 
10+0-05 
F0.05 
0.10 
1 7Q+0.O5 
l-'"-0.05 
-,+0.05 
-0.05 

-2 20+0-05 
^•^"-0.00 
-,+0.05 
0.10 
0.05 



-0.40 
-0.90 



0.10"' 



-0.60^ 



0.20"* 



-1.50] 

-1 90+0-0-"'' 

0.10 
-0.50 
0.05 
0.05 
+0.05 
0.05 
+0.05 
0.00 
0.05 
-0.05 
0.05 
0.00 
0.05 
-0.05 
+0.05 
0.05 
+0.05 
0.05 



-1.70 

-1.10 
-1.60 
-2.20 
-1.70 
-2.20 
-1.40 
0.10 
0.30 



-2.20 



+0.06 
-0.00 



0.044 
0.041 



0.048 
0.004 
0.004 

0.004 
0.004 
0.006 
0.008 



0.005 
0.006 
0.005 
0.004 

0.083 

0.004 

0.058 

0.062 

0.007 
0.012 
0.049 

0.008 
0.009 
0.004 



100 



008 
007 
007 



008 



006 



070 



004 



010 
007 



570 



200 
004 

028 
005 



032 



42 



72 -0, 

-0, 



85 
90 

56 
47 



32 



Table 4 — Continued 



CN Spec 



CN Spec Full Spec UBV^ MSTO'' LIR° HLit'^ BC96« LI^ 



VJK8 LIR lit^ LIR'^ 



HLit<^ 



Object 




Af 


56 




[Fe/H] 


[Fe/H] 


A 


Lge 


Age 


Age 


Age 


Age 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 


[Fe/H] 




Q 


200 


fO 
-0 


100 
050 


^•*J'J-o.io 


r, r,r, + 0.05 

^•^"-0.00 





124 


0.130 


0.631 










-1.00 






NGC2136 





100" 


fO 
-0 


050 
010 


-1 60"^!;-?^ 

^•^^ — 0.10 


-2 20+°-?,^ 
^■^"—0.00 





058 


0.040 


0.081 










-0.40 






NGC2147 





020" 


fl) 

-0 


010 
005 


^•^^-0.05 


-1.80_o 05 





021 




















NGC2156 





100" 


fO 
-0 


050 
010 


-1 60+0-0^ 
^■"'^-0.05 


-1.80_o 05 





077 


0.060 


















NGC2157 





040" 


fO 
-0 


010 

010 


-1 ^o+o-05 

J^-O'J-0.05 


-0 70+0-05 

"■"^-0.05 






0.030 


















NGC2159 





100" 


fO 
-0 


050 
010 


-1.80_o 05 


-i.5oin^ 





056 


0.060 


















NGC2164 





100" 


fO 
-0 


050 
010 


-i.80_o 05 


-o.8ol°;i° 





044 


0.050 


0.058 










-0.60 






NGC2214 





040" 


fl) 

-0 


010 
010 


^•■3U_o 05 


-0 70+0 1° 

'^■"^-0.05 





038 


0.040 


0.052 










-1.20 






SL153 





200" 


fO 
~0 


100 

050 


-2 20+°-^° 
^■^"-0.00 


-i.ooin^ 





054 




















SL492 





010" 


fO 
-0 


005 
001 


3Q+0-05 
'-'•'^'-'-0.05 


-o.solH^ 





005 




















SL538 





006" 


fO 
-0 


001 

001 


W+0-05 
'^•^'^-0.05 


40+0-05 





006 




















SL586 





002" 


fO 
-0 


001 
001 


U-OU_o 05 


r, r,r, + 05 

^■^"-0.00 
























SL601 





004" 


fl) 

-0 


001 
001 


-0 so+°-i° 


_-, -p+0.05 
^■'^-0.05 





004 




















SL639 





0041° 


001 
001 


-1 90+°°^ 
1.9U_o 05 


9 9r,+0.05 





030 


0.015 



















SMC 



-U.500 



+0.05 
0.05 
+0.05 
0.05 
+0.05 
0.05 
+0.05 
0.05 



40+0-05 
U-W„o.o5 
+0.05 
0.00 
+0.05 
0.05 
+0.05 
0.06 



NGC121 
NGC299 
NGC330 
NGC419 



lo.ooo; 



0.008: 



-0.500 
fO.OOl 
'-0.001 

004+° ""l 
U.UU4_o 001 

1 500+°-25° 



-1.70 
-0.20 
-1.50 
-1.20 



-2.20 
0.00 
-1.40 



13.000 

0.012 

0.013 0.007 

0.670 1 



11.900 



0.025 
950 1.200 



70 



90 



19 

82 
70 



M31 



Bol04 


4.0001°:^°° 


-1 20+°°^ 

J^-^'J-0.05 


-0 90+° °^ 
U.MU_o 05 


10.000 








B0II8 


13.000li:°°° 


-1 "S0+°°'5 
i-au_o 05 


40+° °'^ 
>J.4U_o 05 


10.000 








B0II9 


lo.oool-Z 


-i-2ol°:J° 
-0.801°:°;: 


50+° °° 
^■'^^-0.05 


10.000 








B0I24 


10.000lJ:000 


40+°-°'5 
'J-^'^-o.io 


10.000 








Gl 


O.UUU_o 500 


'^•'^'^-0.1)5 


-0 20+°-°'5 
'J-^'J-0.05 










G78 


3.000l°joo 


-0 50"'"° °5 
'^•3'J-0.05 


-n in+° °5 

'^•^'^-0.05 

-n 20+° °5 

'J-^'J-0.05 










G158 


18.000li:0~ 


'J-J^'J-0.05 


10.000 


19.800 


-0.26 




G165 


15.0001?:*° 


-n ^o+o-o'^ 


50+° °° 


10.000 




-0.56 


-0.69 


G168 


20.0001?:°™ 


-0 50"'"°-°5 

'^•'^'^-0.1)5 


50+"-°° 
*J-3'J-0.05 


10.000 


6.200 




-1.13 


G169 


7.oooi;-»» 


, nn + 0.05 
J^-'J'J-0.05 


-1 00+°°5 
^^•'-"^-0.05 


10.000 


7.500 


-1.18 


-1.35 


G170 


20.000lS:°°° 


0.00"^°°^ 


20+° °^ 
'J-^>J_o 05 










G174 


8.0001°:^°° 


S0+° °5 
U.CSU_o 05 


10+°°^ 
U.lU_o 06 


10.000 


7.200 


0.29 


-0.13 


G175 


20.0001°:°°° 


1.0U_o 10 


U.OU_o 06 


10.000 


1.300 




-1.16 



Table 4 — Continued 





CN Spec 


CN Spec 


Full Spec 


UBV^ 


MSTO*" 


LIR<= 


HLifl 


BC96'' 




V.IKS 


LIR lit^ 




HLit^ 


Object 


Age 


[Fe/H] 


[Fe/H] 


Age 


Age 


Age 


Age 


Age 


[Fe/H] 


[Fc/H] 


[Fc/H] 


[Fe/H] 


[Fe/H] 


G177 
G184 
G185 
G190 
G222 


20.000i°:»o 
20.000l«:«S 

4 000+1-°°° 
11 000+°-^°° 


00+° °'5 
-1 "S0+°-°5 
-1 00+°°'^ 
-0 90+° l° 


'-'■^^-0.05 
-0 40+° °'5 

-0 60+° °^ 

-0 60+° °5 

-0 40+° °^ 


10.000 
10.000 
10.000 
10.000 








19.800 
19.800 
6.500 


-0.15 
-1.08 


-0.32 
-0.91 
-1.19 
-1.96 









Note. — All ages are in Gyr. 
''Santos et al. (1995) 
''Hodge (1983) 

•^Thc Lconardi & Rose (2003) comparison literature values come from Olszewski et al. (1991), Cohen (1982), Seggewiss & Richtler (1989), Sagar & Pandey 
(1989), and Piatti et al. (2002). 

Homogeneous lit values (Santos & Piatti 2004) 

<'BC96 fits to colors in BATC filters (Jiang et al. 2003) 

'Line indices (Huchra et al. 1991) 

sCohen & Matthews (1994) 

''Equivalent width of Ca triplet (Olszewski et al. 1991) 



Table 5. Clusters in Red vs. Blue Spectrum Effects Study 



ID 


Loc 


Blue CN 
Spec Age 


Red CN 
Spec Age 


Break CN 
Spec Age 


Entire CN 
Spec Age 


Blue Full 
Spec Age 


Red Full 
Spec Age 


Break Full 
Spec Age 


Entire Full 
Spec Age 


Min Lit 

Age 


Max Lit 
Age 






U.UU4 


on nnn 


U.UU4 


on nnn 
zu.UUU 


A nnn 
4.UUU 


1 f^nn 


fi nnn 
o.UUU 


o nnn 


io.i 




NGC 1783 


LMC 


0.010 


1.500 


9.000 


1.500 


6.000 


1.200 


20.000 


1.200 


0.200 


1.548 


NGC 1831 


LMC 


1.200 


0.600 


1.000 


0.800 


0.400 


0.200 


1.000 


0.400 


0.190 


0.501 


NGC 1847 


LMC 


0.001 


0.040 


0.040 


0.020 


0.080 


0.040 


0.040 


0.040 


0.016 


0.021 


NGC 1856 


LMC 


0.040 


0.600 


0.040 


0.800 


0.600 


1.200 


0.200 


0.400 


0.120 


0.151 


NGC 1866 


LMC 


0.004 


1.200 


0.060 


0.200 


0.100 


0.800 


0.200 


0.600 


0.080 


0.302 


NGC 1868 


LMC 


0.008 


0.400 


2.000 


3.000 


0.200 


1.500 


2.000 


0.800 


0.330 


0.933 


NGC 1978 


LMC 


0.004 


20.000 


10.000 


4.000 


4.000 


20.000 


10.000 


9.000 


2.200 


2.570 


NGC 2004 


LMC 


0.002 


0.020 


0.004 


0.010 


0.020 


0.004 


0.010 


0.010 


0.008 


0.028 


NGC 2070 


LMC 


0.004 


18.000 


10.000 


0.060 


0.008 


0.400 


0.004 


0.006 






NCG 2157 


LMC 


0.001 


1.200 


0.040 


0.040 


0.100 


1.000 


0.060 


0.400 


0.030 


0.030 


NGC 330 


SMC 


0.001 


0.020 


0.020 


0.004 


0.040 


0.004 


0.040 


0.040 


0.007 


0.025 


NGC 419 


SMC 


0.010 


2.000 


3.000 


1.500 


0.200 


3.000 


3.000 


20.000 


0.670 


1.950 



Note. — 



All ages are in Gyr. Fitting ranges are: Blue = A <4000 A, Red = A >4000 A, Break = 3750-4250 A, Entire = full available range. 



